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THEORY  OF  INDUCED  SPATIAL  INCOHERENCE 


I.  INTRODUCTION 

High  gain  laser  fusion  requires  a  highly  symmetric  pellet  implosion,  with  ablation  pressure 
nonuniformities  of  no  more  than  a  few  percent  around  the  spherical  surface  1 2  For  directly  -driven  pel¬ 
lets.  this  symmetry  must  be  achieved  by  the  combination  of  nearly  uniform  laser  illumination  and 
thermal  diffusion  in  the  ablating  plasma.  '  Thermal  diffusion  alone  is  ineffective  in  reducing  lower 
spatial  frequency  nonuniformities  to  acceptable  levels,  especially  at  the  short  optical  wavelengths 
required  for  good  laser-target  coupling  4 ' 

Theoretical  studies  have  shown  that  acceptable  spherical  illumination  uniformity  can  be  achieved 
by  overlapping  a  limited  number  (  2 20)  of  focused  beams,  provided  that  each  individual  beam  profile 
is  smooth  and  controllable. *' g  Earlier  efforts  to  obtain  such  profiles,  however,  have  been  frustrated 
by  the  inherent  imperfections  in  high-power  multistage  laser  systems.  The  cumulative  effect  of 
numerous  small  amplitude  and  phase  aberrations  (both  linear  and  nonlinear)  introduced  by  each  opti¬ 
cal  element  of  a  multistage  laser  produces  large  random  aberrations  in  the  output  beam,  and  hence 
large  random  intensity  nonuiufomnues  at  the  target  surface.  Efforts  to  control  laser  aberrations, 
using  ultra  high  quality  optics  and  extensive  beam  relaying,  have  not  been  completely  successful, 
especially  at  high  energies  and  shorter  laser  wavelengths.10  11  Nonlinear  optical  techniques,  such  as 
phase  conjugation12  and  Raman  beam  cleanup. 13-15  are  effective  at  generating  beams  capable  of  pro¬ 
ducing  nearly  diffraction— limited  focal  spots.  However,  these  techniques  are  not  generally  effective 
in  eliminating  large  residual  intensity  nonuni formities  in  the  quasi  near-field  of  the  laser.1'1'"'  where 
the  pellet  would  have  to  be  placed  in  order  to  obtain  the  required  spot  size  with  a  lens  of  reasonable 
focal  length. 

One  promising  solution  to  the  uniformity  problem  is  the  Induced  Spatial  Incoherence  (ISI)  con¬ 
cept.1 7  ":i  In  the  convenuonal  version  ’  20  of  this  technique,  a  broadband  laser  beam  (bandwidth 
»  1  / 1 puisc)  is  sliced  into  an  array  of  small  beamlets  by  an  orthogonal  pair  of  echelon  structures, 
which  impose  a  different  time  delay  at  each  step.  (The  idea  is  illustrated  in  two  dimensions  in  Fig 
1.)  If  the  delay  increments  to  ■  -  r,  are  chosen  somewhat  longer  than  the  optical  coherence 

time  tc  «  1/Av,  the  beamlets  become  mutually  incoherent.  These  beamlets  are  then  overlapped  onto 
the  target  by  a  lens  of  focal  length  / .  Each  of  them  will  independently  focus  to  the  same  smoothly- 
varying  far-fieid  diffraction  profile  of  total  width  2/X/D,,  provided  that  the  initial  width  O,  is  small 
in  comparison  to  the  aberration  scalelength  r*  in  the  incident  laser  beam.  One  has  the  option  of 
either  centering  all  of  the  profiles  at  a  single  point,  or  (by  slightly  tilting  adjacent  echelon  steps  in 
opposite  directions)  centering  them  equally  around  four  nearby  quadrature  points  in  order  to  control 
the  shape  and  width  of  the  composite  profile.11-20  (Fig.  2)  If  the  incident  beam  aberration  also  con¬ 
tains  some  short  scalelength  components  u  <<  £>,).  which  might  for  example  arise  from  hard 
apertures  or  damage  spots,  the  resulting  energy  will  tend  to  focus  well  outside  the  main  lobe  ot  the 
far-field  diffraction  pattern.  This  component  will  therefore  miss  the  target,  and  will  be  of  little  conse¬ 
quence  as  long  as  it  represents  only  a  small  fraction  of  the  total  energy 
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Although  the  overlapped  beamlets  produce  a  complicated  interference  pattern  at  any  one  instant, 
their  mutual  incoherence  ensures  that  this  pattern  will  evolve  randomly  in  times  of  order  i  It  the 
target  responds  hydrodynanucally .  in  times  ~fKytko  »  it  will  effectively  ignore  this  rapidly  - 
shifting  structure,  and  respond  only  to  the  smooth  tune-averaged  diffraction  profile  For  example,  an 
optical  bandwidth  -  30  cm'1  (easily  achieved  in  Nd  glass  or  KrF  lasers)  provides  /  *  I  ps. 
whereas  time  scales  for  gross  hydrodynamic  motion  of  a  large  high  gain  pellet  shell  are  typically  -  I 
to  10  ns. 

The  idea  of  slicing  the  laser  light  into  beamlets  that  are  subsequently  integrated  at  the  target 
plane  also  forms  the  basis  of  several  other  techniques  for  controlling  irradiation  profiles.  In  the 
random-phase  shift  (RPS)  technique  developed  by  Kato  and  Mima.  22  22  the  beam  slicing  is  performed 
by  a  transparent  phase  mask.24  which  randomly  imposes  a  fixed  phase  shift  of  either  0  or  t  on  each  of 
the  resulting  beamlets.  The  RPS  technique  has  the  advantage  that  it  does  not  require  a  broadband 
laser,  and  the  phase  mask  may.  in  some  existing  laser  facilities,  be  easier  to  utilize  than  ISI  echelons 
However,  the  fixed  phase  relationship  among  the  beamlets  (even  with  broadband  light)  ensures  that 
the  interference  pattern  will  remain  stationary  throughout  the  pulse  duration.  The  lower  spatial  tre- 
quency  components  of  this  pattern  could  not  only  create  nonuruformities  in  the  ablation  pressure,  but 
could  also  seed  self  focusing  or  other  instabilities  in  the  underdense  plasma. 

Several  optical  integrating  techniques  have  been  developed  to  produce  a  "flat-topped"  intensity 
profile  at  the  target.  This  can  be  accomplished  by  an  array  of  prisms. 25  tilted  mirrors.26  2’  or  lens  sys¬ 
tems2*  that  combine  the  beamlets  in  their  quasi-near  fields.  With  these  techniques  one  can  avoid  most 
of  the  deleterious  low  spatial  frequencies  in  the  interference  pattern  if  one  also  uses  low  F'<  number 
optics  (i.e.,  large  angles  between  the  beamlets).  Alternatively,  one  can  achieve  full  ISI  operation  by 
using  a  broadband  laser  and  allowing  the  usual  delay  increments  between  the  tilted  echelon  steps,  as 
illustrated  in  Fig.  3.  All  of  these  quasi  near-fteid  schemes  have  the  drawback  that  Fresnel  diffraction 
of  the  sharp-edged  beamlets  will  introduce  gross  nonuniformities  unless  the  Fresnel  numbers  are 
>  100.  This  requires  very  low  F /number  optics,  which  would  be  unsuitable  for  a  reactor  design."  In 
principle,  both  Fresnel  ripples  and  short  scalelength  intensity  aberrations  can  be  at  least  partially  aver¬ 
aged  out  by  offsetting  the  beamlets  from  one  another  at  the  target.2*  However,  it  is  not  clear  that  this 
averaging  process  could  maintain  adequate  uniformity  at  the  target  if  there  were  any  significant 
amount  of  aberration  in  the  incident  laser  beam.  The  Fresnel  number  requirement  would  be  less 
severe  if  high  power  apodization  were  applied  to  the  edges  of  each  step:  one  possible  technique  for 
achieving  this  is  currently  under  development. 2V  Thus,  the  quasi  near-field  version  of  ISI  may  become 
a  viable  alternative  that  would  be  particularly  useful  for  flat  foil  acceleration  experiments.  The 
remaining  issues  would  then  be  the  limited  depth  of  focus  [approximately  the  spot  size  x  (F/numberi|. 
and  the  cost  and  complexity  of  the  apodized  steps. 


This  paper  presents  a  theoretical  description  of  the  conventional  (far-field)  version  of  ISI;  "flat 
topped"  ISI  will  be  the  subject  of  a  later  paper.  Section  11  describes  the  ISI  configuration  in  greater 
detail,  and  calculates  the  ideal  profiles  in  the  single-focus  and  quadrature  modes  described  above. 
Section  in  examines  nonuniformities  in  the  average  intensity  and  ablation  pressure  due  to  residual 
interference  among  the  beamlets  within  finite  averaging  times.  These  results  are  generalized  to  a  par¬ 
tial  ISI  configuration,  where  beamlets  at  larger  crossing  angles  are  allowed  to  remain  mutually 
coherent  throughout  the  pulse.  Because  the  resulting  interference  pattern  will  contain  only  high  spa¬ 
tial  frequencies,  which  can  presumably  be  smoothed  out  by  thermal  diffusion,  partial  ISI  may  be  one 
possible  way  to  accommodate  a  very  large  number  of  beamlets  without  significantly  affecting  the  tem¬ 
poral  behavior  of  the  pulse  Section  IV  examines  the  perturbing  effects  of  laser  aberration  and  beam- 
let  divergence  on  the  ISI  profile.  Section  V  shows  numencal  simulations  comparing  ISI  results  (aver¬ 
aged  over  times  r  -  100  i,  i  with  those  obtained  by  ordinary  quasi  near-field  illumination  and  bv  the 
RPS  technique.  In  Sec.  VI.  we  present  some  recent  ID  numencal  simulations  indicating  that  (SI 
operation  suppresses  filamentatinn  in  the  underdense  plasma.  Finally.  Sec.  VII  reviews  these  results, 
along  with  those  of  recent  target  interaction  expenments  using  ISI.  and  bneflv  discusses  a  promising 
new  technique  for  implementing  ISI  without  using  echelons.21 
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II.  ISI  CONFIGURATION 

Consider  a  collimated  laser  output  beam  propagating  along  the  + direction.  Its  instantaneous 
field  amplitude  at  point  x  ■  u  .v  )  within  the  transverse  plane  :  =  :L  can  be  written  in  the  form 

Et  (x.n  =  Mix)  F(r)  e"""'.  (2.1) 

where  the  complex  amplitudes  A ( x)  and  Fu)  describe  the  transverse  spatial  structure  (e.g..  due  to 
beam  aberration)  and  the  temporal  dependence,  respectively,  and  b  is  a  constant  to  be  specified  later. 
For  a  broadband  Q -switched  pulse,  where  Av  »  l/t^.Fu)  is  well  approximated  by  a  quasi- 
stationary  stochastic  variable.  Its  correlation  function  Fu  +  n  F*u )  will  exhibit  a  smooth,  localized 
t  -dependence 


FU  +  t)  F*U)  =  I  FU)  | :  7(f) . 


(2.2a) 


where  >(0)  ■  l  and 

7(f)  —  0  for  | f  |  >  ft.  =  l/Av,  (2.2b) 

while  the  mean  square  amplitude  |E(r)|:  describes  the  slow  temporal  behavior  of  the  average 
intensity  over  intervals  -f,*,,*  »  tc.  For  any  quantity  G(t )  =  G(FU)),  the  brackets  GU)  denote  an 
ensemble  average,  which  is  equivalent  to  the  time  average 


e<  Ut' 

GU)r»[  GU')~ 


over  interval  r  in  the  limit  where  rf/r  —  0  and  r  «  The  effects  due  to  averaging  each 

member  of  the  ensemble  over  a  finite  number  of  coherence  times  will  be  considered  in  Sec.  III. 


An  orthogonal  pair  of  echelon  structures,  each  having  Ns  »  l  steps  of  width  £>,.  slices  the 
beam  into  an  Ns  x  .Vs  array  of  square  beamlets.  (If  the  incident  beam  is  circular,  then  only  tA/s:/4  of 
those  beamlets  will  contain  any  appreciable  energy.)  The  nth  beamlet.  which  is  defined  by  the  two- 
dimensional  vector  n  =  (nt.fiv)  with  nr  and  nv  ranging  over  1.2..  ...Vv.  is  centered  at  position 
%  -  (•t„.yB),  where 


",  -  T <*s 


l) 


-  t'Afv 


l) 


(2.3) 


and  delayed  by  time  f„.  The  total  field  amplitude  can  thus  be  written  as  the  two-dimensional  summa¬ 
tion 


Es  (x.f) 


bV  u 
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FU  -  tn)  /t(x)  exp  [i  k0n  ■  (x  -  xj  -  i  uiU 


(2.4a) 


bZ  “ 


C  x  -  x. 


D, 


M, +  (x-xj  (V^), 


+  ~  <X  -  v  ),]f(f  -  t„)  exp  (i*0„  (X  -  x„;  -  (u)f  +  i* „J .  (2.4b) 


where  k  *  2t/X,  utx  i  is  a  “top  hat''  function  equal  to  1  for  t  s  1/2  but  zero  otherwise.  ■  *ir„ 
and  0n  are.  respectively,  the  phase  shifts  and  deflection  angles  imposed  upon  the  beamlets  by  the 
echelon  steps,  and  (  )„  means  that  the  enclosed  quantity  is  to  be  evaluated  at  The  second  order 


expansion  (2.4b)  assumes  that  .-t(x)  vanes  slowly  across  each  beamlet.  with  the  transverse  gradient 
terms  ( V  ,-t  )„  and  <  V  x  v  ^.-n.  descnbing  the  incident  beam  aberration.  For  the  unaberrated  terms, 
we  can  specify  the  normalization 

£  >.4n|!  =  1  (2.5) 

n 

without  any  loss  of  generality.  No  attempt  is  made  (nor  should  be  made)  to  control  the  optical  path 
delays  crn  to  tolerances  of  order  \.  so  the  phase  shifts  are  assumed  to  be  mutually  independent  and 
randomly  distributed  between  0  and  2ir.  The  angular  deflections  0„  are  allowed  up  to  2NS  indepen¬ 
dent  adjustments;  however,  the  most  useful  combinations  appear  to  be  either  0„  =  0  for  all  n  (single¬ 
focus  configuration)  or  0n  alternating  among  the  values  (0.0), (0, -0).( -0.0). ( -d.  -»i)  for  the  quadra¬ 
ture  configuration. 

A  single  lens  of  focal  length  /  now  superimposes  the  beamlets  at  a  target  plane  located  a  short 
distance  ;(«/)  in  front  of  the  focal  point.  The  complex  amplitude  at  this  plane  is  found  by  substi¬ 
tuting  expression  (2.4b)  into  the  diffraction  integral:30 


k  r  * 

£(x.z.r)  =  - — - -  \d-x'  Es(x'.t)ex p  / 

2ir(f  -  z)  J 


.  it  lx*  -  *  I  ' 


Uf  -  z) 


.  it  lx' I - 


(2.6a) 


k  b  exp  fi  x(x,z ,/ )]  „  ,  ....  *xn  ‘  x  ! 

- JZJ-7) £  FU  -  ij  exp  JZ)  -  «  —  , 


lond^  Id/ 2d*y  '4»  +  :  « 


x  exp  X  -  J  x„  -  (/  -  Z)  0„  ■  S  +  i  rfj r  y  1 

where  x  (x.z.t)  *  it(/  -  ;)  -wt  +  k\  x|  :/2 (/  -z).  and 


<t>' „(c )  m  *„  +  kz  I  x„ |  :/2 /  if  -z) 

is  a  random  phase  factor.  In  writing  expression  (2.6).  we  have  ignored  the  propagation  distances  /„ 
between  the  echelon  steps  and  the  focusing  lens.  This  is  justified  as  long  as  the  target  lies  in  the 
quasi  far-field  of  the  beamlets;  i.e.,  |  z  |  must  remain  small  in  comparison  to  the  beamlet  focal  depth: 

| z  |  <<2X/:/Df  (2.7) 

Expression  (2.6)  then  takes  advantage  of  the  fact  that  the  far-field  amplitudes  are  independent  of 
except  for  small  phase  curvature  terms  kln\x\2/2fz.}0  The  resulting  phase  differences 
k(ln  “  /,  )  I  x|2/2/2  between  beamlets  are  negligible,  so  these  terms  can  be  removed  from  the  summa¬ 
tions  and  incorporated  into  x(x.  z .  t )•  Because  |$|  £  D,/2,/2.  the  |£|:  term  in  the  exponential  of 
(2.6b)  can  be  ignored  as  long  as  condition  (2.7)  is  satisfied.  One  then  obtains 


Dxb  expfi x(x.c./)|  „  kx„  x 

E(x.z.i)  =  - T.  FU  -  f„)  exp  <*„(:)  -  /— - 


X 


+ 


(?^>™'  VX 


(T  j.  V  X/1  )„  :  717i  sine  x 


*  -  ( z/f )  *„-(/-■ 


where 


sine  (0) 


sin  9X  sin  6y 


a  a  [f  —  z )  X/ D  |  —  f\/D\ 


(2.10) 


are,  respectively,  the  functional  form  and  effective  width  of  the  far-field  diffraction  pattern  due  to  a 
single  unaberrated  beamlet  of  width  D , .  The  condition  required  for  negligible  separation  of  the  beam- 
lets  due  to  their  angular  divergence  is  |x„|  \z  \  /fa  «  1,  or  in  terms  of  the  F/number  / /D  of  the 
total  aperture  D  =  NSDU 


f  ID  >>  z/2a  . 


(2.11) 


The  largest  defocusing  1 2  |  that  would  be  required  in  a  spherical  illumination  geometry  would  be 
the  pellet  radius  ;  hence,  condition  (2.11)  can  be  adequately  satisfied  by  using  high  F/number 
optics.  Finally,  we  note  that  the  far-field  condition  (2.7)  can  be  rewritten  as  //D,  >  >  \z  i  /2a.  and 
hence  follows  immediately  from  (2.11)  because  D/D,  =  Ns  >  1. 

Returning  to  Eq.  (2.8),  we  expand  the  sine  function  up  to  second  order  in  x„z/f  a.  approximate 
/-2s*/,  and  finally  specify  b  *  (8r/cDr)l/2: 


E(x.z.t)  = 


exp(i*(x.  z,  t)l  £  F<{  ~  r«)  exP  ' 


kx„  -  x 


x  <  A 


”  ^  /  fc”  *"  (  V  )"  •  7  A  V  X  r  S'nC  »  - - 


V  >n  +•  IT  I'  W^n 


(2.12) 


Expression  (2.12)  contains  all  of  the  information  needed  to  generate  not  only  the  ideal  IS1  profiles, 
but  also  the  deviations  arising  from  multi-beamlet  interference,  aberration,  and  beamlet  divergence 
effects. 

To  complete  this  section,  we  evaluate  the  ideal  ISI  profiles  for  the  single-focus  (i.e.  complete 
overlap)  and  quadrature  configurations  described  previously.  ISI  generally  requires  all  of  the  beam- 
lets  to  be  mutually  incoherent;  i.e.,  the  correlation  function  defined  in  (2.2)  must  satisfy 


7(f„  -  f„.)  =  5„  „• 


(2.13) 


In  the  scheme  proposed  here,  this  condition  is  achieved  by  choosing  the  incremental  delays  to 
satisfy 


4 


'•;?£i 


At  >  tc  =  \/Av 


(2.14) 


for  all  n.  [A  factor  At/tc  =  1.5  will  usually  '«  sufficient  to  ensure  that  (2.13)  is  well  satisfied  if 
Av  =  Av(FWHM).]  The  ideal  profiles  are  then  found  from  expression  (2. 12)  by  applying  (2.13)  to  the 
ensemble-averaged  intensity  /  (x.z.t))  *  (c/8t)<  \E(x.zj)\2 \  and  discarding  the  aberration  and 
beamlet  divergence  terms.  For  the  complete  overlap  configuration  (all  =  0)  one  obtains 


where 


f,0'(x.t)  =  (l/a:)  S(irx/a) 


P*(t)  *  E  \F(‘  ~  r-)l:  M.l: 


(2.15) 


(2.16) 


is  the  ensemble-average  of  the  total  power,  and 


S{0)  *  sinc:(0)  = 


sin  9,  sin  9r 


(2.17) 


is  the  diffraction  profile  of  a  single  beamlet.  Approximately  82%  of  the  incident  energy  is  contained 
in  the  central  lobe  of  S(xx/a).  as  defined  by  -a  s  t.  y  s  a.  The  quadrature  configuration 
corresponds  to  x  ■  0„  =  +/3(-d)  for  even  (odd)  values  of  nt  and  y  0„  =  (-d)  for  even  (odd) 

values  of  /t, .  This  gives 

PQ\x,t)  =  P6’(t)  S  |  rx  ~  &  V  — —  1  (2.18) 

a'Z  L  a  a  J 

where  6a  =»  fl3  is  the  focal  spot  offset,  a  and  »  both  assume  the  values  ±1.  and  the  quantities 


Po*(t)  «  £  £  |F(r  ~  («)|-  M, 


/>,f-(r)  -  £  £  IF(r  -  f,)l:  !  etc.  (2.19) 


refer  to  summations  over  the  even  (e)  and  odd  (as  steps.  The  optimum  value  of  6a /a  depends  on  the 
specific  target  application.  For  the  broadest  possible  flat-topped  profile,  as  shown  in  Fig.  2b.  the 
optimum  choice  would  be  6a  =  0.43a ;  to  optimize  the  illumination  uniformity  and  coupling  in  a 
multi-beam  spherical  geometry,18 19  the  best  choice  centers  around  oa  ~  0.22a . 

In  order  to  avoid  significant  temporal  broadening  effects  due  to  the  echelons,  the  total  differen¬ 
tial  delay  should  remain  short  in  comparison  to  the  puisewidth:  i.e. . 


(mu  (min  33  At  <<  tg 


(2.20) 


[The  combination  of  this  criterion  and  condition  (2.14)  thus  places  restrictions  on  the  maximum 
number  of  steps  and  the  minimum  optical  bandwidth.]  If  (2.20)  is  applicable,  then  IFu  -rB)i: 
may  be  removed  from  the  summations  in  (2.16)  and  (2.19),  which  are  then  subject  to  the  normaliza¬ 
tion  condition  (2.5);  e.g.,  the  total  average  power  can  be  written  as 

Po(t )  =  I  F(r )  I :  (2.21) 

In  the  quadrature  configuration,  the  summations  will  satisfy  the  condition  /V*u>  +  /Y~<n 
+  Pq*U)i  +  Pq-(t)  =  P0(t)  ;  thus,  in  the  ideal  case  where  these  summations  can  be  exactly  bal¬ 
anced,  the  profile  becomes 


/'e'(x,r)  =  (l/a:)  Pt)(t)  Q  ( xx/a  ,6a  /a  ) . 


(2.22) 


■- 


where 


.  1  _  _  _  t  -  ilou  \  -  i 'ha 

Qixx/a  .oa  /a)  =  —  X'  x-  5  x - .  x - 

4  “*  —  a  u 

M  ✓ 

The  5< t  x/a )  and  Q(rx/a  .oa  /a)  profiles  are  compared  in  Fig.  2. 


III.  NONUNIFORMITIES  DUE  TO  RESIDUAL  INTERFERENCE 


The  target  responds  to  quantities,  such  as  intensity  or  ablation  pressure,  averaged  over  some  fin¬ 
ite  number  r/t  of  coherence  times.  Unlike  ensemble  averages,  these  finite  time  averages  will  retain 
small  nonuniformities  due  to  residual  interference  among  the  beamlets.  even  if  the  beamlets  are  mutu 
ally  incoherent.  This  section  will  examine  those  nonuniformities  in  the  limit  where  perturbations  due 
to  aberration  and  beamlet  divergence  are  negligible.  For  simplicity,  we  consider  only  the  complete 
overlap  configuration;  the  quadrature  configuration  will  produce  similar  results  for  the  cases  of 
interest  where  5a  <  0.4a.  The  amplitude  Eix.t)  can  then  be  approximated  bv  the  zeroth  order  term  in 
Eq.  (2.12),  with  0n  =  0.  Thus,  the  time-averaged  intensity  profile  can  be  written  as  the  multimode 
expansion 


/(x.r ) .  ■  <c  /8t)  E  (x.r ) ; ;  . 


1  X 

-  —  s  (t— )  v'  Pm  (t)  ,  exp  (iK*,  X) . 

a  m 

where  S  ( rx/a )  is  the  ideal  diffraction  envelope  defined  in  Eq.  (2. 17).  and 


PmV)  r  ■  E  AaA*„.m  exp  -  «<*'„  F(t'  -  tn)  F*U'  ~  .  (3.2a) 


K„,  m  k(x„  -  x„)/f  =  m  kDjf  =  (2T/a)m  *  m  AAT  (3.2bi 

are.  respectively,  the  complex  amplinides  and  spatial  frequencies  of  the  transverse  modes  The  m  =  0 
amplitude  Pn(t )  r  =  F(rM: .  is  the  time  average  of  the  total  power,  so  1 3 .2a  1  represents  a  straightfor¬ 
ward  generalization  of  (2.16).  All  of  the  m  *  0  modes,  which  describe  the  interference  pattern,  are 
statistically  independent,  except  for  the  constraint  P-m(n  .  =  Pm u)  *  required  to  ensure  that  /  1  x.m  . 
is  real. 

The  actual  target  acceleration  is  determined  by  the  average  ablation  pressure  pix.n  .  .  rather 
than  the  intensity  itself.  In  general,  p(x.n  T  can  be  related  to  the  instantaneous  intensity  /ix.n  by  a 
space-time  convolution  over  the  conduction  zone  between  the  ablation  surface  and  absorption  region. 
This  relationship  will  depend  upon  the  detailed  hydrodynamic  behavior  within  that  zone,  and  upon  the 
optical  wavelength,  average  irradiance  level,  and  other  parameters  that  also  affect  the  absorption  in 
the  underdense  plasma. 

To  a  good  approximation,  one  can  ignore  the  temporal  convolution  and  relate  the  time-averaged 
pressure  p(x.t)  ,  to  a  spatial  convolution  over  the  time-averaged  intensity  /< x.n  .  in  the  limit  where  r 
satisfies  the  condition 

^  pulse  f hydro  —  r  f  *Fh  -  (3.3) 

where  r,  and  t,h  are,  respectively,  the  ion-acoustic  transit  and  thermal  diffusion  times  across  the  con¬ 
duction  zone.  These  times  are  given  by  the  approximate  expressions 


where  d  is  the  effective  width  of  the  conduction  zone,  c,  =  9.8  x  103  (Z/.4)1'2  [7,(eV)]l/2  (with  atomic 
number  .4  and  charge  state  Z)  is  the  ion-acoustic  velocity  in  cm/sec,  vf  =4.2  x  I07  [7,(eV)|l/2  is  the 
electron  thermaJ  speed.  r„  =3  4  x  I03  [7, (eV))5'2  mf  In  A  is  the  electron-ion  collision  time31  in  sec,  nr 
is  the  electron  number  density,  and  InA  is  the  Coulomb  logarithm.  For  recent  NRL  planar-target 
experiments  using  527  nm.  2  ns  pulses  at  intensities  - 1013  -  I014  W/cm1.  20  we  estimate 
T,  -500  eV  d  -50  »m.‘  and  n,  -  1022  cm"';  these  parameters  then  yield  r,  -320  ps  and  t,h  -360  ps. 
Condition  i3 .3)  is  thus  reasonably  well-satistled  as  long  as  the  gross  hydrodynamic  response  of  the 
target  occurs  on  nanosecond  time  scales.  A  -.tmilar  result  is  expected  for  reactor-sized  targets  operat¬ 
ing  at  shorter  optical  wavelengths,  longer  pulsewidths.  and  higher  intensities.  For  example,  for  a 
possible  high  gain  pellet  driven  by  a  few  MJ.  1.4  ^m  laser  with  a  5-10  ns  pulse  at  5  x  1014  W  cm2,  we 
estimate  7,  -  ’50  d  -  20  -  50  nm .  and  n,  -  4x  I022  cm-5:  this  gives  r,  -  100-250  ps  and 
r,„  -  80  -  500  ps. 

The  spatial  convolution  over  /(x.n  .  is  approximated  here  by  using  the  heuristic  "cloudy  day 
model1 2  to  account  for  thermal  smoothing  effects  within  the  conduction  zone.  If  p  iK./i  .  and 
T  iK.r  ,  are  the  transverse  spatial  Fourier  transforms  of  pix.n  .  and  Itx.n  ..  respectively,  then  the 
cioudv  day  model  provides  the  simple  result 

p  ‘  K.i  I  ,  =  C  /  i  K..r )  .  exp  <  -  K  d  \ .  1 3  4 1 

where  C  is  a  constant.  Taking  the  Founer  transform  of  Fq.  (3.1)  and  recalling  definition  i2  I’i.  one 
obtains 

p  i  K./ 1  .  *  C  ^  SiK  -  K,^)  Pm  1 t )  .  exp  -  K  d  i .  1 3  5 1 


where  C  *  C /a:.  and 

5  iK)  ■  fd2r  5  it  —i  exp  <  -i  K  \) 

J  a 

1 1  -  A,  /A A  Ml-  A.  AA  hi  for  A,  and  A  <  A  A 
~  0  for  A,  or  A,  >  A  A 

is  a  triangular  function  whose  effective  width  AA  =»  lv/a  =  kD,/f  is  the  intermode  spacing.  |See 
Eq.  (3.2b).] 

We  are  primarily  interested  in  focal  spot  diameters  <  -  a  \  that  are  large  in  comparison  to  the 
width  of  the  conduction  zone,  specifically  the  case  where 

AAd  =  2t  d la  <  I  (3.6) 

Expression  (3.5)  is  therefore  well-approximated  by 

P  <K.n  ,  »  CT5(K  -  K*,)  Pm  in  T  exp  <-  1  K„.d). 

m 

which  transforms  back  to 

pix.n  .  =  C  "  S(t  —  i  v  pm  in  .  exp  (  -  Km  I  d)  exp  n  K,  -  xi.  (3.7) 


Comparing  this  result  to  expression  (3.1).  we  see  that  thermal  smoothing  tends  to  filter  out  the  higher 
spatial  frequency  modes  arising  from  the  interference,  but  has  little  effect  on  the  diffraction  envelope 

In  the  narrowband  limit  where  tc  exceeds  the  total  time  delay  -  tmin  =  ,VS:  _u ,  the  beamlets 
would  remain  mutually  coherent.  This  limit  corresponds  to  the  RPS  technique." :3  The  amplitudes 
Pm(t))T  retain  large  randomly-phased  contributions  for  all  m  *0,  in  addition  to  the  m  =  0  term 
P0(t)  T  =  ,  |  F  (r )  |  -  r;  hence  the  average  intensity  profile  (3.1)  would  consist  of  a  stationary  but  ran¬ 
dom  multimode  interference  pattern  modulated  by  the  smooth  S(xx/a)  envelope,  as  illustrated  in  Fig. 
1.  Similar  considerations  apply  to  the  ablation  pressure  profile  (3.7),  except  that  the  higher  spatial 
frequency  components  of  the  interference  pattern  are  suppressed  by  the  exp  (-|K„|d)  factors.  This 
filtering  effect  thus  provides  the  primary  basis  for  the  RPS  technique.  In  the  following  paragraphs, 
we  will  derive  approximate  expressions  for  the  RMS  mode  amplitude  of  the  nonuniform  time- 
averaged  intensity  and  ablation  pressure.  These  expressions  can  be  used  to  evaluate  the  standard 
deviations  of  the  interference  nonuniformities  for  both  the  narrowband  and  broadband  cases. 

A.  Complete  ISI  Operation 

For  complete  ISI  operation,  Ai>  must  be  broad  enough  to  satisfy  conditions  (2.14)  and  (2.20). 
thus  ensuring  mutual  incoherence  for  all  of  the  beamlets.  The  m  *  0  components  of  Jm(t)  T  will 
then  decrease  in  relation  to  70(/)>r  as  tc/r  decreases;  i.e.,  expressions  (3.1)  and  (3.7)  will  become 
dominated  by  the  S  (rx/a)  profile  as  the  interference  structure  fades  away.  In  the  ensemble  limit 
where  tcir  —  0,  expressions  (3.2a),  (2.13),  and  (2.5)  yield 

Pn(t)r-  PmU)  r  «  <W  |f(0|'  r.  (3.8) 

and  both  ;7(x,/))r  and  p  (x.t)  r  (as  given  by  the  approximation  (3.7)]  reduce  to  the  exact  5(irx/a)  pro¬ 
file.  The  double  average  in  (3.8)  allows  an  averaging  over  the  slow  nonrandom  variation  of  P0(i) 
in  cases  where  r  could  become  a  non-negligible  fraction  of  f^. 

The  relative  magnitude  of  the  interference  components  of  / (x.r )  ,  and  p  (x.t)  r  can  be  quanti¬ 
fied  by  comparing  the  RMS  mode  amplitudes  for  m  *  0  to  that  of  the  DC(m  =  0)  component. 
Expressions  (3.1)  and  (3.7)  suggest  the  most  appropriate  definitions  are  (for  m  *  0) 

a,(m)  rn  \  pm(t),r\l  [r- /  P0(t)  *  1/2  (3.9a) 

ap(m)  m  \  Pm(t ))T  | 2  1/2  exp  ( -  j  K*,  j  d)/  /*,)(/) >rl/2  =  <7,<m)exp<  -  |  m  i  AKd)  (3.9b) 

where  the  ' )  brackets  again  denote  an  ensemble  average,  and  [recalling  (3.2b)]  AK  =  k  Dx/f  -  2 r/a. 
The  total  RMS  nonuniformities  due  to  all  of  the  statistically  independent  m  *  0  modes  are  then 
evaluated  from  the  vector  sums 


"1/2 

2  <jf(m)  I 

"  ,  1  1/2 

oP  *  |  ^  <V(m> 

(3. 10a. b) 

m  *  0 

L  j 

m  »  0  ' 

L  J 

A  subsequent  paragraph  will  show  that  P^t)<;  1/2  =  P, thus,  expressions  (3  10a. b)  may  be 
regarded  to  a  good  approximation  as  the  standard  deviations  of  the  spatial  nonuniformities. 


From  expression  (3.2).  we  obtain 


!  P»d)  t  I  =  Z  Z  '*•  '*%*«  4%  -  m  exp  [/(o'.  -  o'...  -  o',  -t-  o'.  .m;j 

•  m' 

X  f'"'  —  f'*'  —FU'-IJ  F*(t'  -tn^.m)F*{t"  -f.)  FU”-ta^m)  (3.11) 

Jr  r  r 

for  ail  values  of  m,  including  m  =  0.  We  now  assume  that  the  broadband  light  is  sufficiently  chaotic 
to  satisfy  Gaussian  statistics,  thus  allowing  the  factorization32 

Fit' -tjF*{f  -  F*  U"  -i„)F  (t"  - 1„  rm) 

-  F(»'-i.)P(r'-f,.J  F*U"  -t„)  F(t" 


+  F  (t'  ~  tj  F*  d"  -  fn ) 


(3.12) 


Temporal  mode  locking  behavior  due  to  “instantaneous"  nonlinear  effects,  such  as  self-focusing  or 
harmonic  conversion,  will  invariably  have  some  effect  on  the  statistics.  However,  the  Gaussian 
model  is  expected  to  remain  a  reasonable  approximation  as  long  as  the  number  of  temporal  modes 
'puIk/'c  (  =  remains  large  and  the  average  nonlinear  phase  shift  does  not  exceed  2t  by  a  large 

amount.  Substituting  Eq.  (3.12)  into  (3.11)  and  recalling  expression  (2.2a)  and  condition  (2.20),  one 
obtains 

.  r  -*•  r  w#  •  .  i  +  r  •  • 

|F(r')|-’'[  ~  IFO”)!2' 

x  £  £  +  *  «•«  exp  [l  (0'.  —  0'.*.  —  <t>\.  +  <j>'n.  )] 

■  a' 

X  [7ft..*  ~  fj7*(f«'+„  -  In)  +  7ft'  ~  t"  -  t.  +  f.)  7*(f '  -  I"  -  t..m  +  /.  ..)]  (3.13) 
Equation  (3.13)  gives  for  the  m  =0  term  [recall  the  normalization  (2.5)1 

dt '  i  „  .  i  r'  ~T  dt " 


Pod)  l  =  l  —  \F(t')\ 2  |  —  [ F (/") | : 

Jl  T  T 

1  +  Z  L  1  A*\z  ^»'l:  I  yd'-t"  -  f.  +  /„•) I  - 1  . 


(3.14) 


The  second  term  within  the  square  brackets  is  nonvanishing  only  when  t'  -  t"  -  r.  +•  r..  »  0;  hence, 
if  condition  (2.20)  is  satisfied,  one  can  approximate  |  F(t')  | :  by  |  F(t")\ 2  when  evaluating  the  con¬ 
tribution  due  to  this  term.  In  the  limit  where  the  total  time  delay  /Vs2  At  satisfies  additional  criterion 
Afs2  A t  <  <r,  condition  t'  -  t"  -  r.  +  r,.  ■  0  can  be  attained  for  all  n  and  n'  over  most  of  the 
(t.t  +  r)  interval.  The  time  integrals  will  then  be  approximately  independent  of  n  and  n',  so  the  sum¬ 
mations  can  be  evaluated  immediately  with  the  aid  of  Eq.  (2.5).  Thus. 

Pod) )r  *  |F(f)l:  ;  +  \F(I)\1  2  ,  7iMS  tc/T, 

where  |F(r)|2  varies  negligibly  within  times  <  <  r,  and 
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f  - 

»  t 


y<t ) 1 


(3  15) 


10 


is  a  constant.  In  the  opposite  limit  where  r  <  <  (Vs2  A/ .  the  second  term  of  (3. 14)  becomes  compar¬ 
able  to 


i  F(/) )  -  -  r  V  \An  \A  -  \FU)\-'-r/Ni . 

n 

which  is  again  small  as  long  as  Nf  >  >  1.  Recalling  expression  (3.8),  one  thus  obtains  the  expected 
result 


P0(t)  2  a  |  Fit)  |2  2  =  P0u)  r  2 

in  all  cases.  The  standard  deviation  in  the  case  where  r  >  A/s;  a/  is  approximately 


(3.16) 


<J0  m 


PoU))T- 


1/2 

/  P o(0  r 


S  7RMS(ft  /r)l/2  <  1/Afj  • 


(3.17) 


assuming  that  r  is  sufficiently  short  in  comparison  to  rr.„  that  the  slowly-varying  quantity  | Fit)  l: 
satisfies  i(  |F(r)  |2'2  rl/2  =  |  Fit)  |2  In  contrast  to  Eqs.  (3.9)  and  (3.10),  this  quantity  measures 

temporal  fluctuations  of  the  entire  profile  amplitude,  rather  than  spatial  nonuniformities. 


For  m  *  0,  expression  (2.13)  eliminates  the  7(fn*m~tii)  terms  in  Eq.  (3.13),  and  the  random 
phase  factors  <t>'„  -  <t>'a  +  m  -  <*>'„•  +  +  m  ensure  that  only  those  contributions  with  n'  =  n  will  add 

coherently  in  the  second  set  of  terms.  Equation  (3.13)  then  reduces  to  the  result 


l^-(O),  I2:  =((  I F  (/)  | 2  )2  )T  7gMS  —  ^  I  I2  M, 


(3.18) 


Substituting  expressions  (3.18)  and  (3.16)  into  (3.9)  and  (3.10),  while  again  assuming  tc  «r« 
one  obtains  the  general  results  for  m  *  0 

<7,(m)/a,  (m)  *  a, /a,  *  ap  (m)/ap(m)  =  ap/ap  *  7rms  (rc/T)l/2,  (3.19) 

where  the  factors 


or(m)  m 

£  1  |  2  |  /4n  *  m  |  2 

1/2 

a,  m 

£  ff,2(m) 

n 

m  *  0  j 

ap  (m)  m  af(m) exp  (-  |  m  |  &Kd) . 


£  ^(™) 
m  *0 


1/2 


(3.20a.b) 


(3.21a,b) 


describe  the  nonuniformities  in  the  narrowband  (RPS)  limit.  Expressions  (3.19)  and  (3.21)  show 
explicitly  that  two  independent  mechanisms  are  combining  to  smooth  the  ablation  pressure.  The  ther¬ 
mal  diffusion  factors  ap(m)  and  ap  filter  out  the  higher  spatial  frequencies,  but  have  little  effect  on  the 
longer  wavelength  modes  unless  A  Kd  -  l.  The  temporal  factor  7rms0c/t)i/2  *  (rc /r)l/2,  which  also 
smooths  the  irradiation,  reduces  the  relative  magnitude  of  the  nonuniformities  as  the  number  of 
independent  random  intensity  contributions  r/tc  increases;  it  has  the  exact  form  that  one  would  expect 
from  a  random-walk  process  where  the  RMS  nonuniformities  increase  as  ( r/tc)u 2,  while  the  average 
fluence  increases  as  r/rc. 

Useful  approximations  to  (3.20)  and  (3.21)  can  be  obtained  using  the  normalization  condition 
(2.5),  provided  that  the  laser  output  beam  is  not  too  badly  aberrated  (e.g.,  peak/vailey  spatial  intensity 
variations  <  2:1).  In  the  case  of  a  nominally  flat-topped  circular  beam  of  diameter  D=/VSDI.  one 


can  aproxinuue  |4„|:  by  the  average  value  4/tjVs:  for  |  x„  |  sD/2,  and  zero  elsewhere.  For 
.V5  »l,  the  n  summation  in  (3.20a)  is  therefore  well-approximated  by  <4/ ir/ Vf):5,  where  S  is  the 
area  of  the  overlap  region  between  two  circles  of  diameter  Ns  whose  centers  are  separated  by  distance 
I  m  |  *  0;  thus,  one  obtains 

a,(m)  »  (8,/:/tJVs)  t/2  -  sin'1  ( |  m  |  /N,)  -  (  |  m  j  //Vs)(  1  -  |  m  |  :/Nf)'/l  ,  (3.22) 

which  is  plotted  in  Fig.  4. 

The  standard  deviation  of  the  integrated  intensity  can  be  found  immediately  by  combining  (3.19) 
and  (3.20a,b)  and  using  the  normalization  (2.5): 

°i  •  (1  -  E  Mml4)'/2  -  (1  -  4/x  N’)ul  <  1  (3.23a) 


a!  <  7RMS  (tc/T) 


1/2 


(3.23b) 


Expression  (3.23a)  shows  that  without  temporal  averaging,  the  RMS  amplitude  of  the  irradiance 
nonuniformities  remains  nearly  independent  of  the  number  of  steps  Ns,  and  is  comparable  to  the  DC 
term.  By  increasing  Ns,  one  merely  shifts  the  spectrum  to  higher  spatial  frequencies,  without  signifi¬ 
cantly  affecting  the  total  energy  in  the  interference  pattern. 

The  standard  deviation  of  the  ablation  pressure  can  be  approximated  by  substituting  expression 
(3.22)  into  Eqs.  (3.21a,b)  and  (3.19);  thus  ap  -  yrms  (ic/t)uz  a„,  where  the  2D  summation 
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r  N$ 


L 

*  *  0 


T 

y  -  sin 


1/2 


j*l  (1  -  exp(-2 AKdNsItl)  i  .  (M  -  m /Ns  s  1) 


(3.24) 


is  plotted  vs  AKdNs  in  Fig.  5.  To  illustrate  these  results,  we  consider  the  parameters  used  for  the  ISI 
simulations  shown  in  Sec.  V;  i.e.,  X  =  527  nm,  /  =  6  m,  Ns  =  16  steps  across  the  FWHM  of  the 
beam,  stepwidth  D|  =  1  cm,  and  absorption-ablation  distance  d  =  50  ^m.  One  then  finds 
A Kd  *  {k  D\/f)d  *  1  and  &KdNs  =  16,  thus  giving  op  =  6.6%  from  Fig.  5  and  op  <  0.66%  for 
r  >  100  cc . 

For  short  wavelengths  and  high  gain  targets,  where  condition  (3.6)  is  well-satisfied  (e.g., 
AKd  -  India  s  0.5),  Fig.  5  shows  that  op  depends  essentially  on  the  single  variable  AKdNs ,  and  fol¬ 
lows  the  simple  approximation  op  <  2u2/HKdNs .  Thus,  for  a  wide  range  of  conditions  applicable  to 
high  gain  laser  fusion,  the  standard  deviation  of  the  ablation  pressure  can  be  approximated  by 


Op  <  7 RMS  (tc/T)U1  Op  *  O,  at 


(3.25a) 


<  2,/VA KdNs  -  a  /2  /z  rdNs  =  (X/2,/z  xd)  ■  (f /D) 


(3.25b) 


where  f /D  -  /  /Ns  D,  is  the  overall  F/number  of  the  beam.  As  a  numerical  example,  we  choose 
Ns  m  60  steps/beam,  a  ■  1.5  mm.  d  =  20  to  50  ^m,  coherence  time  tc  =  1  ps.  and  averaging  time 
r  ■  1  ns.  (The  1.5  mm  would  allow  pellets  up  to  2  mm  radius  with  quadrature  illumination,  and  the 
20  to  50  ntn  is  an  appropriate  range  of  absorption-ablation  distances  for  1/4  nm  light.)  The  standard 
deviations  for  the  narrowband  and  broadband  cases  are  then  ap  < 


11  to  28%  and  op  <  0.4  to  0.9%, 


respectively,  and  the  total  F/number  is  /  /Ns  D ,  =*  a/Ns  X  =  100  for  X  =  0.25  ^m. 


We  close  this  discussion  with  a  caveat  concerning  the  slow  residual  fluctuations  in  the  ideal 
ablation  pressure 

p0(x.i))T  a  C'  P0(t))T  S(vx/a)  (3.26) 

due  to  the  P0(t )>T  fluctuations.  According  to  the  discussion  leading  to  expression  (3.17),  the  standard 
deviation  <r0  of  P0(t))T  is  the  lesser  of  (rr/r)1'2  or  \/Ns,  and  can  therefore  amount  to  several  percent, 
even  for  r  -  1  nsec.  This  poses  no  problem  as  long  as  one  uses  only  a  single  laser  beam.  In  a  multi¬ 
beam  spherical  illumination  geometry, however,  it  could  lead  to  a  small  illumination  imbalance, 
unless  all  of  the  beams  exhibit  nearly  the  same  random  temporal  behavior.  One  can  avoid  the  prob¬ 
lem  by  using  a  single  master  oscillator  and  balancing  the  net  time  delays  within  the  different  chains  as 
closely  as  possible  (at  least  to  within  times  «  t,h).  This  will  allow  some  interference  to  occur 

among  beamiets  arising  from  different  beams;  however,  the  large  angles  (typically  >  20°)  between 
these  beams  insure  that  the  resulting  transverse  wavelengths  will  be  on  the  order  of  microns,  and 
hence  too  short  to  cause  any  problem. 

B.  Partial  ISI  Operation 

Complete  ISI  operation  may  be  difficult  to  implement  in  a  large  aperture  or  highly  aberrated 
laser  because  the  large  number  of  incoherent  beamiets  that  could  be  required  would  invalidate  condi¬ 
tion  (2.20).  The  number  of  beamiets  is  determined  by  the  step  width  D ,  ,  which  depends  upon  the 
target  diameter  <2 a  -  2X//D,  and  the  amount  of  aberration  in  the  laser.  For  X  <  1/2  ^m,  reason¬ 
able  focal  lengths  (  e.g.  /  5  30  m),  and  high  gain  targets  of  a  few  mm  diameter,  one  would  typically 
need  D,  <  1  cm.  With  a  100  cm  aperture  per  beam  this  would  require  /Vx-  >  104  beamiets;  it  would 
significantly  affect  the  temporal  behavior  of  the  pulse,  because  for  tc  -  l  ps.  Ns:  A/  >  Nf  tc  -  10  ns. 

One  possible  way  to  avoid  this  problem  is  a  partial  ISI  scheme,  in  which  the  echelon  step 
sequence  would  be  repeated  after  every  N  <  Ns  steps  across  the  aperture,  (e.g..  see  Fig.  12). 
Expression  (2.13)  would  then  be  replaced  by  the  more  general  condition 

7<'l.  ~  '«•)  “  £  <Vn'  -  Vi..  (3.27) 

m 

where  (mt,mv)  range  over  all  positive  and  negative  integers,  including  zero.  Because  the  repeated 
steps  (m  *  0)  remain  coherent,  they  create  stationary  interference  patterns  at  the  target.  However,  if 
one  chooses  S  sufficiently  large,  these  patterns  will  have  only  high  spatial  frequencies  of  magnitude 
N&K.  2 ,/ZN±K,  2 N&K,  etc,  and  can  therefore  be  smoothed  out  by  thermal  diffusion.  The  relative 
magnitude  of  the  interference  components  of  p  (x,/»r  can  again  be  calculated  from  (3.9)  and  (3.10). 
using  expression  (3.27)  and  arguments  similar  to  those  presented  before.  Thus.  (3.19)  and  (3.21) 
generalize  to 

'  ,  tc  ,  ,  1 

op  s  7 ksts—  T,  °p(m)  +  E  op-(JVm)  1  (3.28) 

T  m  »  0  m  *  0 

where  the  n  summations  again  extend  over  all  ,VS;  beamiets.  The  first  term  within  the  brackets 
describes  the  residual  incoherent  contributions,  and  is  identical  to  expression  (3.19)  and  (3.2lbi;  the 
second  describes  the  contribution  of  stationary  interference  patterns  due  to  the  repeating  steps.  If  one 
chose  N  =  1  (corresponding  to  the  random  phase  case)  then  the  second  term  would  be  identical  to  the 
thermal  smoothing  factor  (3.21b),  and  would  become  dominant  for  tc  «  t.  If  v  is  chosen  large 
enough  to  satisfy  the  condition  2N<lKd  »  I.  then  the  second  term  becomes  negligible  and  (3  28) 
reduces  to  the  usual  ISI  result.  Thus,  for  a  large  high  gain  pellet  where  condition  (3  6)  is  well- 
satisfied.  the  standard  deviation  of  the  ablation  pressure  would  still  be  described  by  (3  25a. bi  The 
ISI  configuration  recently  proposed  for  the  Sirius-M  test  reactor  design1"  provides  an  instructive 


numerical  example.  Here,  each  1.2  m  aperture  (square)  beam  was  divided  into  a  total  of  ,V5  =  240 
steps,  which  repeat  after  every  60  steps.  With  a  stepwidth  D,  =  5  mm.  focal  length  /  =  30  m.  and 
wavelength  \  =  1/4  jim,  this  gives  a  -  f  X/D,  =  1.5  mm.  Assuming  (as  before)  that  d  -  20  to  50 
ntn  and  r/tc  -  1000,  we  calculate  ar  -  0.1  to  0.25%  from  (3.25a.b).  The  overall  F/number  of  each 
1.2  m  beam  is  F/25. 

IV.  ABERRATION  AND  BEAMLET  DIVERGENCE  EFFECTS 

This  section  examines  perturbations  of  the  ideal  target  profile  due  to  laser  aberration  and  beam- 
let  divergence.  In  order  to  separate  out  these  effects,  the  analysis  will  be  carried  out  in  the  ensemble- 
average  limit  tc/T  —  0,  assuming  complete  ISI  operation  defined  by  conditions  (2.14)  and  (2.20). 
Because  the  resuiting  expressions  contain  no  interference  terms,  and  therefore  involve  only  the  long 
scalelength  functions  such  as  S(irx/a)  and  its  first  and  second-order  gradients,  they  should  apply 
equally  well  to  the  average  intensity  or  the  ablation  pressure  as  long  as  condition  (3.6)  is  satisfied. 

A.  Completely  Overlapped  Beamlets 

The  complete  overlap  configuration  is  obtained  by  setting  all  /3„  —  0  in  Eq.  (2.12).  Applying 
conditions  (2.13)  and  (2.21)  to  (2.12).  and  retaining  only  perturbation  terms  up  to  second  order,  one 
finds  the  average  intensity  l(x.z)  =  (c/8ir>t  E(x.z)  I :  : 

I(x.z)  -  ~  £  f  M.l2  ~  Re(/tn*  B„)  ■  +  i  B„Bb*  :  V,  j  5  j  t- 

a‘  ■  a  2  j  d 


sine 


[Re  04*  V 


(  ?  x/t  )„(  V  ^,4*)J  :  V  x  V  sine  j  ir  — 

a 


(4.1) 


where 

B(x)  ■  (/ Hknv  ^.4)  -I-  (; //) x.4  .  (4  2) 

[Recall  that  (  )„  evaluates  the  enclosed  quantity  at  point  x„.]  The  slow  temporal  variation  of 

/  (x.c.n  and  P^n  (on  a  scale  of  t and  the  c -dependence  of  B„<;)  have  been  suppressed  here  in 
order  to  simplify  the  notation.  Because  ,4<x»  varies  slowly  across  the  beamlet  width  £>,.  the  2D  sum¬ 
mations  over  n  can  be  approximated  by  integrals.  For  example  the  normalization  (2.5)  can  be 
replaced  by 

*  J  ,4(x>i;-^-r  =  l,  (4  3) 

n  Df 

and  integration  by  parts  can  be  used  to  eliminate  the  last  pair  of  terms  in  (4  1);  i.e.. 

£  Re</1*  V  V  .4  1,  »  -7  f  — 7  (A*  V  V  A  +C.C.) 

•  -  u  \ 

j2  x 

-  -  f  — t  iV  A)  (  V  (^)  a  -  TiV  Am  V  A*) 

:  O;  A  7 

Equation  (4  1)  therefore  reduces  to 

Itx.Z)  =  [1  -  C ,  i  c  )  V_V_|  /""(X).  (4  4) 

where  t',h<x i  =  i  Pt,u: i  Sitx.ui  is  the  ideal  target  profile  defined  in  Sec  II.  and 


(4.5) 


C,u:)  ■  C,U)-f^-BB* 

J  Df  ’  J  Df 

are  the  aberration  coefficients.  These  coefficients  can  be  rewritten  in  terms  of  the  phase  <t>(x)  and 
modulus  |/4(x)|  of  the  amplitude  A  (x);  i.e.. 


C,(z)  =  |/t(x)l:[/0(x)  +  ^x], 

J  Of  / 

h(z)  =  M(x)|:f[/tf(x)  +  ^x|2  +  [^Vx/«  | a  |]2!. 


(4.6a) 


(4.6b) 


where  tf(x)  a  ( 1  /A: )  v  X4>(x)  is  the  refraction  angle  acquired  by  the  aberrated  wavefront  around  point  x, 
and  the  squared  vectors  [  ]2  in  (4.6b)  represent  diadic  products  [  ]  [  J.  Expressions  (4.6a.b)  can 

be  interpreted  physically  as  intensity— weighted  averages  of  the  bracketed  quantities  over  the  laser 
aperture. 

It  is  useful  to  relate  C,(z)  and  C2(z)  to  spatial  moments33  of  the  coherent  (but  aberrated)  focal 
distribution  /c(x,z))  that  the  laser  would  produce  without  ISI.  As  shown  in  the  appendix,  expression 
(4.6a)  is  just  the  centroid  of  the  coherent  beam: 

C,(z)  =  |xc).  a  |  Ic(x.z))xd2x /f(Ic(x,;)  d2x  (4.7) 

This  result  can  be  readily  understood  in  terms  of  geometrical  optics.  A  ray  passing  through  point  x  in 
the  laser  aperture  will  intercept  the  target  plane  at  point  xc(x)  =  fd(x)  +  xz  // ;  expression  (4.6a)  then 
averages  all  such  contributions  over  the  entire  laser  aperture.  For  the  second  moment  |xcxc|,  which 
is  a  measure  of  the  RMS  width  of  ;/c(x,z)),  geometrical  optics  would  average  the  diadic 
+  xz//l2  over  the  laser  aperture,  thus  giving  the  first  set  of  terms  in  (4.6b).  A  more  thorough 
treatment  including  diffraction  yields  the  complete  expression  (Appendix) 

C2(z)  =  |XCXC|.  a  j  fc(x,z)  xxd2 T/J  rc(x.z)>dz. x.  (4.8) 

although  the  geometrical  optics  contribution  will  usually  remain  dominant  in  saturated  lasers,  where 
most  of  the  aberration  resides  in  the  phase. 

The  aberrated  ISI  profile  (4.4)  can  now  be  expressed  in  terms  of  the  coherent  moments: 

/(X,z»  »  [1  -|Xcl.-  +  jtXcXci.  :  VAV  J  /,0,(x)  (4.9) 

By  dividing  out  the  maximum  value  /<03(0)  >,  one  can  rewrite  this  in  a  dimensionless  form  suitable  for 
comparing  the  magnitudes  of  the  perturbation  terms:  i.e., 

/(x,z))/'/<0’( 0)  »  S(rx/a)  -  t(|xc|.  /a  )■  S'frx/a  ) 


where 


+■  —  t2  (!xc*clf  /«') :  S”(irx/a). 


S'la)  a  V(I5(a),  S"  a  V„V„S (a). 


(4.10) 


(4.11) 


■av.Vv;, 


and  ■  d/da.  The  functions  S(at.O),  S',(at,0),  and  .S"„(at,0)  are  compared  in  Fig.  6.  Equation 
(4.9)  can  also  be  expressed  in  an  alternative  form  by  expanding  /(0’(x)  around  point  |xc|.  and,  as 
usual,  retaining  terms  only  up  to  second  order  in  the  gradients;  i.e., 


Hx.z)>  =  [1  +  -^-|AxcAxc|.  :  V  x  V  J7<0,(x  -  |xc|. ). ,  (4.12) 

where 

|AxcAxc|.  =  Ixcxc|;  -  ixc|.  lxc|; .  (4.13) 

Because  /,0l(x)>  is  symmetric  around  x  =  0,  it  is  apparent  that  |xc|  also  defines  the  centroid  (and 
hence  the  alignment)  of  the  ISI  profile.  A  nonzero  value  of  |xc).  would  be  of  little  consequence  in  a 
flat  foil  experiment,  where  precise  alignment  is  unnecessary.  However,  in  a  spherical  geometry 
where  each  laser  beam  must  be  properly  centered  onto  the  pellet,  |xc|.  should  be  treated  as  a  pertur¬ 
bation  term  [c.f.  Eq.  (4.10)]  if  it  cannot  be  entirely  removed  by  the  alignment  system. 

If  the  aberrations  are  spatially  random,  one  can  use  a  simple  statistical  model34  to  estimate  both 
the  magnitudes  and  reproducibility  of  the  perturbation  amplitudes  in  Eq.  (4.10).  This  model  should 
provide  a  reasonable  description  of  perturbations  arising  from  small-scale  effects  such  as  turbulence, 
multiple  surface  nonuniformities,  and  residual  lens  aberrations,  which  contribute  the  lion’s  share  of 
the  far-fieid  broadening  and  are  the  most  difficult  to  eliminate.  Systematic  lens  aberrations,  such  as 
astigmatism  in  angularly— multiplexed  KrF  systems,3536  can  be  largely  compensated  within  the  laser. 
Large-scale  effects,  such  as  thermal  gradients,  tend  to  cause  gradual  beam  steering  effects  rather  than 
any  significant  broadening,  and  can  be  balanced  out  in  real-time  by  an  automatic  alignment  system. 
For  simplicity,  the  model  will  treat  only  phase  aberration,  which  is  the  primary  culprit  responsible  for 
both  beam  steering  and  broadening  effects  in  real  lasers.  The  random  phase  aberration  is  character¬ 
ized  by  an  effective  coherence  width  <  D  (where  D  is  the  laser  aperture),  over  which  the  RMS 
phase  shift  is  *.  The  local  refraction  angle  9<x)  -  (1  /k)V  ^4>(x)  can  then  be  characterized  using  a 
number 


=  (D/sa>)z  >  1  (4.14) 

of  statisticaily-independent  contributions  9,  whose  ensemble-average  values  are  9  =  0  (assuming  that 

the  laser  is  properly  aligned)  and  whose  RMS  amplitudes  are 


10,  I  2  =  T/fcS^  =  \/2s„b  . 


7KMS  ”  !  »/  I  “  h 

From  Eqs.  (4.6a,b)  (4.7)  and  (4.8),  one  thus  obtains  for  the  relevant  moments 


(4.15) 


r2 


Ixci 


f  '  * 

=  2  L  9i  ■  l*C*cl  =  ~T  L9I9I 


N. 


ab  j  »  I 


^  ah  f  *  I 


where  the  subscript  ;  has  been  dropped  in  order  to  simplify  the  notation. 


(4. 16a.b) 


Combining  Eqs.  (4. 14)-(4. 16a)  and  using  the  statistical  conditions 
9j  ■  9j  ■  0jw«6(/,  one  obtains  the  ensemble-average  |xci;  =  0  and  the  RMS  centroid 


=  0  and 


■  l*cl!'  '  “  -  f9RMS ab  -  f\/2D  :  (4  17) 

i.e..  the  RMS  misalignment  due  to  small  scale  beam  steering  effects  is  approximately  half  the 
diffraction-limited  far-field  width  f\/D.  (This  expression  probably  represents  an  upper  bound  on  the 
misalignment  because  at  least  some  of  the  contributions  to  it  may  be  correctable  by  an  automatic 


alignment  system.)  The  RMS  value  of  the  S'(tx/u)  coefficient  in  Eq.  (4.10)  thus  reduces  to  the  sim¬ 
ple  result 


Ri  m  jr(!xc|)ff,M5/ci  =  iir/2)f\/Da  (4.18a) 

=  r/2Ms  ■  (4.18b) 

where  Eq.  (2.10)  and  the  relationship  ,VS  =  D/D,  was  used  to  obtain  (4.18b).  [Note  that  the  criterion 
specified  in  Sec.  II  requires  ,VS/.V(UJ  =  sah/Dt>>  1  (e.g.  >  3)  to  ensure  the  accuracy  of  the  perturba¬ 
tion  expansion  (4. 10).  For  the  two  large-scale  ISI  examples  discussed  in  Sec.  Ill,  where  we  chose 
Nj  =  60  and  240,  the  respective  values  of  /?,  are  2.5%  and  0.6%.  It  should  be  noted  that  expression 
(4.18a)  is  also  applicable  to  conventional  illumination  schemes,6-8  where  a  would  represent  the  diam¬ 
eter  of  the  quasi  near-field  profile  at  the  target.  Obviously,  aberration-induced  beam  steering  is  an 
issue  that  can  affect  all  spherical  illumination  schemes,  not  just  ISI. 

The  mean  value  of  the  second  moment  !xcxc|)  will  reduce  to  its  diagonal  elements 
|.t<f f  =  (vc=|;  =  ||xc|:|)/2  as  long  as  the  aberration  is  spatially-random.  Combining  Eqs.  (4.14), 
(4.15)  and  (4.16b),  one  obtains  for  the  mean  square  diameter  of  the  (non-ISI)  focal  spot 

4t|  |  Xq  I  *1)  ~  (2/^*ms)"  =  (/X/s^)-  “  {f\/D)~N^, ,  (4.19a) 

which  can  be  rewritten  in  the  intuitively  satisfying  form 

N, *  =  2:|  |  xc  |  :|}i/2/(/\/D) 


=  (RMS  diameter  )  /  (  diffraction  -  limited  diameter  ).  (4. 19b) 

The  effective  magnitude  of  the  S "(xx/a)  coefficient  in  Eq.  (4.10)  is  then 

R2  m  (tt/8)4(|  |  xc  |  :l)/u: 


=  (ir/SKiV^Di/Dr  =  (ir/8W,i/A/s:.  (4.20) 

which  tends  to  broaden  the  ISI  profile,  as  shown  in  Fig.  7.  If  the  laser  aberration  were  only  a 
few  times  diffraction  limit,  one  could  virtually  eliminate  this  broadening  term  without  using  an  exces¬ 
sive  number  of  steps;  e.g.,  a  choice  of  iVj/ZV^  >  10  would  reduce  R2  to  about  1%  or  less.  However, 
the  high  power  multistage  lasers  envisioned  as  fusion  drivers  are  likely  to  be  heavily  aberrated,  with 
typical  values  N &  -  10  to  20  for  each  30  cm  of  aperture.  It  is  therefore  expected  that  will  lie 

in  the  range  of  3  to  5,  giving  the  respective  values  R2  =  15%  to  5%.  The  function  S  ’lxx/a)  is 
smooth  and  symmetric  (Fig.  6),  so  a  perturbation  amplitude  of  this  size  will  be  tolerable  as  long  as  it 
remains  reproducible  to  within  -  10%.  Because  of  the  spatial  averaging  process  inherent  in  expres¬ 
sions  (4.6a,b)  one  can  expect  a  high  degree  of  reproducibility. 

The  magnitude  <\R2  of  the  non-reproducible  part  of  the  perturbation  amplitude  can  be  estimated 
from  the  RMS  deviation 


A(|xc|:!  ■  ('.( I  xc  ]  2|:  -  1  |  Xc  1  :l):],/2  (4.21a) 


where  expression  (4.21b)  follows  from  (4.16b)  and  the  statistical  independence  of  the  9,  contribu¬ 
tions.  Under  chaotic  conditions,  one  can  approximate  |0;|4'  -  1 0;  I : :  by  ,  thereby  obtaining 
A|  |  xc  |  :j  =  fl9kus/^ab-  The  resulting  contribution  to  the  perturbation  amplitude  is  then 


A R2  *  (ir/8)4A|  |  xc  |  :|/a: 


=  R'jHan  =  (71^/8  )Njh/Nsz.  (4.22) 

according  to  (4.14),  (4.19a)  and  (4.20).  For  the  simulations  discussed  in  Sec.  V  {Nah  =  5  x  diffrac¬ 
tion  limit  and  Vs  -  16  steps  across  the  FWHM  of  the  beam)  one  expects  A Rz  -  2.4%;  however,  a R: 
will  be  much  smaller  under  the  conditions  applicable  to  a  fusion  driver.  For  example,  a  20  x 
diffraction-limited  beam  in  a  30  cm  aperture  with  60  steps  would  have  A Rz  =  0.7%;  a  1.2  m  com¬ 
pound  aperture  consisting  of  a  4  x  4  array  of  such  beams  (i.e. .  Ntlb  =  80  and  Ns  =240)  would  allow 
only  A Ri  ~  0.17%. 

B.  Quadrature  Configuration 

The  above  perturbation  treatment  gives  similar  results  for  the  effects  of  phase  aberration  and 
beamlet  divergence  on  the  quadrature  profile,  at  least  for  the  cases  of  interest  where  the  offset  oa 
satisfies  5a  <  0.4 a.  Unlike  the  complete  overlap  case,  however,  the  quadrature  profile  /‘°Vx)  may 
be  affected  by  laser  intensity  nonuniformities  even  if  the  usual  perturbation  terms  in  Eq.  (4.10)  were 
negligible.  To  examine  this  effect  in  greater  detail,  we  rewrite  expression  (2.18),  assume  condition 
(2.20),  and  recall  the  normalization  +  P,f  +  p0~*  +  P(f”  =  P «/: 

IiQ\x ))  =  -^f-  Q(xx/a,5a/a)  +  £  £  PS  —  fS - 

^  u.,  a'.**  8 a 


x 


(4.23) 


where  Q(irx/a.5a/a)  is  the  ideal  quadrature  defined  in  Eq.  (2.23).  We  now  expand  the  S  functions 
up  to  second  order  around  S(vx/a)  (valid  for  small  da  /a)  and  combine  the  P$'  .  P,J  ’  summations 
to  finally  obtain 


/'y,(x)i  = 


Q{vx/a  .5 a  /a) 


One  can  null  out  most  of  the  contribution  to  coefficients  (4.25a-c)  just  by  carefully  aligning  and 
orienting  the  echelons  within  the  beam;  however,  in  practical  lasers,  where  the  saturated  near-tield 
profile  varies  during  the  pulse  duration  (and  may  not  even  be  completely  reproducible  from  shot-to- 
shot),  the  residual  contributions  should  probably  be  treated  as  random  variables.  Each  of  them  would 
then  have  the  approximate  RMS  value 

Ckms  =  (— Air5'-)l/‘(A|/t  |  ")rms  *  (1/2^'A/jHA  |  A  |  _)rMS/(  |  A  j  ,  (4.26) 

where  ( M  |  :)av,  *  1/Mf  and  (AM  |  :)RMS  respectively  describe  the  average  and  RMS  deviation  of  the 
intensity  across  the  laser  profile.  For  example,  a  ratio  (AM  |:)RMs/(M  |:)a„  =  0.15  (corresponding 
to  *  2:1  peak7valley  variation)  would  give  CRMS  =  0A/Ns,  which  gives  CRMS  =  0.7%  for  ,VS  =  16 
and  0.2%  for  <V5  =  60. 

V.  NUMERICAL  SIMULATIONS  OF  ISI 

The  beam  smoothing  behavior  of  ISI  can  be  clearly  illustrated  by  numerical  simulations  of  the 
time-averaged  intensity  and  ablation  pressure  profiles  for  a  planar  target.  These  simulations  were  car¬ 
ried  out  for  0.527  /xm  light  by  a  3D  numerical  code  (CHAOS),  which  evaluates  the  diffraction 
integral  (2.6a)  for  an  incident  laser  beam  of  arbitrary  aberration,  but  ignores  refraction  and  filamenta- 
tion  in  the  underdense  plasma.  Thermal  smoothing  effects  are  modelled  by  the  heuristic  “cloudy 
day”  approximation,  as  defined  in  expression  (3.4),  assuming  in  these  examples  an  effective 
absorption-ablation  distance  of  d  =  50  ^m.  For  the  simulations  shown  here,  the  chaotic  temporal 
behavior  of  the  laser  light  is  modelled  by  Gaussian  statistics.  The  real  and  imaginary  parts  of  the 
complex  amplitude  F(t)  are  assigned  independent  random  values  with  a  Gaussian  probability  distribu¬ 
tion,  and  these  values  fluctuate  randomly  from  one  coherence  time  to  the  next.  Some  of  the  simula¬ 
tions  were  also  carried  out  using  other  statistical  distributions,  such  as  random  phase  with  constant 
intensity;  the  results  indicated  that  the  ISI  smoothing  should  be  relatively  insensitive  to  the  statistics  of 
the  incident  light. 

Figure  8  illustrates  the  problem  encountered  by  the  conventional  illumination  scheme. 7IU0n 
where  the  target  is  simply  placed  in  the  quasi  near-field  of  the  focusing  lens.  In  order  to  model  a 
saturated  laser  beam,  the  incident  aberration  was  chosen  to  reside  primarily  in  the  phase:  i.e..  the 
-20  cm  beam  had  a  5  x  D.L.  angular  spectrum,  but  only  -30%  peak/valley  intensity  nonuniformi¬ 
ties.  This  was  accomplished  by  choosing  the  random  complex  amplitude  /l(x>  to  satisfy  Gaussian 
statistics,  then  attenuating  the  intensity  nonuniformities  by  replacing  the  modulus  M  (x>  by 
M(x)|1/8.  The  overall  laser  profile  was  modelled  by  multiplying  the  resulting  values  of  ,-t (x>  by  a 
hypergaussian  envelope  (exp [ — (2  |  x  [  /D)10]j,/2.  At  the  target,  where  the  beam  has  focused  down  to 
S  0.5  mm,  the  irradiance  nonuniformities  are  significantly  worse  because  of  the  phase  aberration  in 
the  incident  beam.1011  Large  random  nonuniformities  are  also  very  evident  in  the  ablation  pressure 
profile,  although  thermal  smoothing  has  removed  the  higher  spatial  frequency  components. 

Figure  9  models  the  original  version  of  the  RPS  technique.”  23  A  pair  or  reflecting  echelons 
(Ns  =  20,  D,  =  1  cm)  have  been  inserted  into  the  beam  with  their  alternating  steps  tilted  to  produce  a 
quadrature  profile  at  the  target  plane,  which  is  now  located  in  the  far  field  of  a  6  meter  lens.  (This 
configuration  is  similar  to  the  one  used  in  recent  NRL  experiments  with  broadband  0.527  ^m  light.30) 
In  this  case,  however,  the  bandwidth  Av  is  assumed  to  be  sufficiently  narrow  that  all  400  beamlets 
remain  mutually  coherent,  thus  resulting  in  a  complicated  interference  pattern  that  remains  stationary 
throughout  the  pulse.  The  thermal  smoothing  mechanism  was  unable  to  filter  out  the  lower  spatial 
frequency  components  in  the  corresponding  ablation  pressure.  The  maximum  peakyvalley  deviation 
from  the  ensemble-averaged  profile  p(x),  was  found  to  be  37%.  which  is  consistent  with  the  RMS 
value  ap  -  6.6%  calculated  in  Sec.  III.  Because  the  RPS  technique  allows  a  smaller  stepwidth  than 
ISI,  it  would  be  possible  to  reduce  this  nonuniformity  by  increasing  the  number  of  steps.  [E.g. .  see 
Eq.  (3.25b).]  However,  it  would  also  be  necessary  to  either  increase  the  laser  aperture  or 
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decrease  the  focal  length  /  in  order  to  maintain  the  same  spot  size  a .  For  example,  a  tenfold 
decrease  in  bp  would  require  Ns  =  200,  which  would  in  turn  require  F/3  optics. 

Figures  (10)  and  (11)  show  the  time-average  intensity  and  ablation  pressure  profiles  for  the  com¬ 
plete  ISI  case,  in  which  Av  is  now  chosen  large  enough  to  ensure  mutually-incoherent  beamlets. 
Although  the  averaging  interval  r  =  100  tc  was  chosen  relatively  short  in  this  simulation,  the  quality 
of  the  target  profiles  already  shows  dramatic  improvement  over  previous  results.  The  residual  low 
amplitude  structure  on  the  average  intensity  profile  has  been  all  but  completely  eliminated  from  the 
ablation  pressure  p(x)  T.  The  maximum  (peakyvalley)  deviation  from  the  ideal  (ensemble-averaged) 
pressure  p(x)  was  found  to  be  2.4%,  which  is  consistent  with  the  RMS  value  <x„  =  0.66%  calculated 
in  Sec.  III.  On  the  basis  of  Eq.  (3.25a),  we  expect  the  maximum  deviation  to  be  <1%  with  averag¬ 
ing  times  r  >  500f<.. 

Figures  12  and  13  simulate  the  partial  ISI  case,  where  the  time  delay  sequence  now  repeats  after 
the  10th  step;  thus  steps  1  and  11,2  and  12.  etc  remain  coherent.  The  resulting  increase  in  fine  scale 
interference  structure  (of  wavelength  \f  / 10 £>,  =32  pun)  is  clearly  evident  in  the  average  intensity,  and 
would  persist  throughout  the  entire  pulse.  As  a  result  of  thermal  smoothing,  however,  this  structure 
is  effectively  removed  from  the  ablation  pressure;  e  g.,  the  corresponding  contribution  to  expression 
(3.28)  is  proportional  to  exp  (  -AvNd/a)  =*  exp  (  -20).  The  3.9%  peakyvalley  deviation  found  here 
is  somewhat  larger  than  in  the  previous  case,  but  that  difference  appears  to  be  statistical. 

VI.  FILAMENT ATION  SIMULATIONS 

The  filamentation  instability  is  an  important  consideration  in  the  use  of  ISI.  Because  the  hydro- 
dynamic  response  time  for  filamentation  (typically  -  100  ps)  is  much  longer  than  the  laser  coherence 
time  (about  1  ps),  it  would  appear  that  filamentation  would  be  easily  suppressed  by  ISI.  In  addition, 
random  density  fluctuations  driven  by  the  constantly  changing  intensity  structure  would  act  to  diffuse 
any  residual  intense  spots,  and  spoil  filamentation  further. 

On  the  other  hand,  the  smoothing  effect  of  ISI  is  incomplete  on  the  filamentation-hydrodvnamic 
time  scales;  time-averaged  intensity  fluctuations  of  -  10%  persist  over  this  period,  and  may  still 
drive  filamentation.  Also,  there  will  be  short  moments  during  which  the  local  intensity  incident  on 
target  will  be  much  higher  than  average.  [For  Gaussian  statistics,  the  intensity  distribution  is 
(l//  )exp <-///  )•  |  Because  filamentation  is  inherently  a  nonlinear  optical  instability,  one  may  argue 
that  these  momentary  intensity  peaks,  rather  than  the  time  averaged  fluctuations,  dominate  the  fila¬ 
mentation  behavior.  Thus  the  effect  of  ISI  on  light  filamentation  is  not  readily  discernible  from  sim¬ 
ple  arguments. 

A  numerical  code  SELFOCT  has  been  developed  to  simulate  thermal  and  ponderomotive  self 
focusing  in  laser-plasma  interactions,  and  particularly  any  filamentation  of  the  light  inhomogeneities 
that  may  occur  under  ISI  conditions.  The  code  is  time  dependent  and  Cartesian  in  two  spatial  dimen¬ 
sions.  Light  propagation  is  handled  by  the  parabolic  wave  equation,  while  the  plasma  dynamics  are 
calculated  by  semi-lineanzed  fluid  equations  combined  with  flux-limited  heat  flow. 
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Laser  propagation  in  SELFOCT  is  described  by  the  paraxial  WKB  approximation  to  Maxwell's 
wave  equation 


(c*  V -  -  dl/dr)  E  =  ( 1  -  iv„  /ui)uj,;£  .  (6  I) 

This  equation  assumes  a  quasi-neutral  plasma  characterized  by  the  local  plasma  frequency 
(j )p  =  u)(n,(x,c.n/ur]l/:  and  electron-ion  collision  frequency  v„  <x{n,ix  ,z.n/Tr'  :ix.z.r  >■''  One  can 
separate  the  field  variations  into  fast  and  siow  space-time  scales  by  writing 
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(6.2) 


Eix  .  z .  t )  =  ♦(*  .Z.t  )exp  [i  j  k0(z  ‘)dz’  -  i  wi  ] . 

o 

where  ♦(x.c.n  is  the  slowly-varying  complex  amplitude,  fc0  =  € ,1/ 2 ( c ) / c  ,  and  =  1  - n,„<; )/nr  is 
the  reai  part  of  the  dielectric  constant  determined  by  the  unperturbed  steady-state  density  neo(z>.  In 
the  paraxial  WKB  approximation,  ♦u.c.n  and  ka(z)  remain  slowly-varying  over  spatial  scales 
\~2x/k0(z).  Substituting  (6.2)  into  (6.1)  and  retaining  only  the  first  order  variations  of  +(.t .z.t)  in  ; 
and  r ,  one  obtains 


j  d  id.  ,  d: 

°  j^dz  +  v„  dt  J  +  ‘  dz  dar 
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=  ~  (€r  ~  M  +<€,)♦• 

C" 


(6.3) 


Here.  er(.r.c,r)  =  1  -  nrl.x  .z.r)/nc  and  €,(x.z.r)  =  v„tx  .c.r)«,(.r  ,z.t)/n[U>  are.  respectively,  the 

real  and  imaginary  parts  of  the  instantaneous  plasma  dielectric  constant,  and  =  c  €,!, - )  is  the 
group  velocity.  This  equation  can  be  simplified  by  applying  the  additional  transformation 


*  =  * 


MO) 
k0(z ) 
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(6.-1) 


where  €«(:)  =  v„(z)neo(z)/nc<Jj,  and  rewriting  all  quantities  in  terms  of  the  time  t  =  / -c/v..u ). 
Redefining  the  axial  coordinate  as  drj  =  dz/i^l(z)  and  normalizing  both  spatial  coordinates  to  the  opt¬ 
ical  wavelength  X  =  2xc  /u,  one  finally  arrives  at  the  parabolic  equation  in  canonical  form 


(4  xid/dy  +  d2/dx2)Mx  .z.t)  -  -4ir  |^6(x.c.t) 
where  €■€,+/€,  is  the  complex  dielectric  constant. 
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Mx  .z.t). 


(6.5) 


Equation  (6.5)  is  solved  using  a  split-step  fast  Fourier  transform  (FFT)  technique.'7  To  accom¬ 
modate  the  FFT  technique,  the  computational  mesh  must  be  periodic  in  t .  and  this  affects  the  manner 
in  which  the  incident  (c  =  0)  ISI  laser  field  is  numerically  constructed.  The  ISI  echelons  separate  the 
broadband  laser  beam  into  a  multitude  of  mutually-incoherent  beamlets.  which  are  then  overlapped  at 
the  target  by  the  focusing  lens.  In  the  SELFOCT  code,  these  overlapping  beamlets  are  approximated 
by  incoherent  plane  waves  propagating  at  slightly  different  angles  0„  =  .t „//,  where  x„  is  given  by 
expression  (2.3)  and  |0„  |  <  <1.  The  incident  optical  field  is  constructed  on  a  calculational  mesh  in 
Fourier  k-space,  with  each  plane  wave  assigned  a  transverse  spatial  frequency 
km  =  MO)  tan  0,  »  ko(Q)0„ .  where  MO)  *  u/c .  Modes  are  exactly  resolvable  when  |  kx„  |  Xle„  is  an 
integer  multiple  of  2r.  where  X,„  is  the  transverse  width  of  the  mesh;  modes  that  lie  between  the 
points  on  this  k-space  mesh  are  assigned  to  the  nearest  mesh  point.  Any  modes  outside  of  the  resolv¬ 
able  frequency  range  (i.e.,  IM  !  X,„  or)  are  assigned  to  the  kr  =0  (DC)  mode.  The  complex  ampli¬ 
tude  of  each  beamlet  is  assumed  to  satisfy  Gaussian  statistics,  as  in  Sec.  V,  with  independent  random 
fluctuations  from  one  coherence  time  to  the  next. 


The  plasma  response  to  the  laser  field  is  found  by  solving  linearized  fluid  equations  that  include 
the  ponderomotive  force  term.  Combining  the  continuity  and  momentum  equations  for  the  fluid 
plasma,  linearizing  in  the  flow  velocity,  and  ignoring  the  second-order  term  V  In (n<  i V <.y  one  obtains 
a  driven  ion-acoustic  wave  equation  for  the  quantity  L(x.z.t)  m  ln(ne): 
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dr  dt 
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where  a  ■  Ze:/4mrm,<jjj,  c,:  =  ( ZTt  +  T,)/m . ,  and  vta  is  a  heuristic  damping  rate.  Linearization  is  jus¬ 
tified  by  the  fact  that  the  fractional  density  and  electron  temperature  perturbations  were  found  to  be 


typically 


especially  in  the  ISI  runs  where  rapid  fluctuations  of  the  interference  structure 


tended  to  suppress  large  perturbations.  In  any  case,  most  of  the  non-ISI  simulations  were  carried  out 
using  a  steady-state  treatment,  which  does  not  require  linearization.  [The  use  of  the  variable  la(nr> 
instead  of  nt ,  preserves  the  correct  ponderomotive  nonlinear  steady  state  behavior.)  The  ion  tempera¬ 
ture  is  assumed  to  be  constant;  ion-electron  energy  equilibration  times  are  typically  on  the  order  of 
nanoseconds  in  these  plasmas. 

The  ion-acoustic  propagation  and  heat  flow  are  expected  to  occur  primarily  in  the  transverse  i.r > 
direction.  This  stems  from  the  following  two  considerations:  (i)  under  the  paraxial  wave  conditions 
assumed  here  the  intensity  variations  occur  primarily  along  .t ,  and  (ii)  because  the  optical  propagation 
speed  ( -  c)  is  so  much  larger  than  cs,  the  ion-acoustic  '•wake"  will  propagate  nearly  normal  to  the 
optical  axis  ; .  Assuming  then  that  all  density  disturbances  propagate  only  along  c ,  with  relatively 
small  variations  in  the  sound  speed,  one  can  Fourier-transform  Eq.  (6.6)  to  yield: 


- r  +  qk  — —  +•  k  ~  Lik  .z.t  )  = 

dt'2  df  I 


1  a  I  *  j  -  +  c2 


where  t'  ■  tc, /X.  k'  m  kt\  and  q  *  v,u/k,i\  is  the  ratio  of  the  imaginary /real  frequency  of  the  ion- 
acoustic  wave.  (In  this  paper  q  is  taken  to  be  1/2).  The  first  term  on  the  RHS  of  (6.7)  is  due  to  the 
ponderomotive  force,  and  the  second  term  is  the  plasma  pressure  due  to  thermal  gradients.  These 
source  terms  are  assumed  to  be  approximately  constant  over  the  time  step  used  in  the  code  (typically 
a  picosecond  or  so).  The  solution  to  (6.7)  for  a  source  constant  from  r0‘  to  r0'  +  r,’  is: 

L(k t0'  t')  -  exp  ( -gk  't ,72)  <{  \/K[dHk '.  r0')/dt()'  +■  qk  'L  ik  ' G ]  sin  (*/,') 

r  ^ 

+  +  C  !  cos  j  -  G  (6.8) 

where  k  m  1 1  -  k '  and  G  ■  la  i  ♦  <£  r0’  *  <  s’/—)  ■ '  '.  r(1'  •+■  f  |72)|,'c,;, .  This  result  is 

inverse  Fourier-transformed  to  yield  the  plasma  density  n,u.c.r„'  -r ,')  at  the  next  time  step. 


acoustic  wave. 


The  electron  temperature  is  found  by  solving  the  relevant  energy  balance  equation. 

3  d  .  H 

2  '  dt  Tf  dx 


where  Qe  is  the  electron  thermal  heat  flux.  H  is  the  joule  heating  source  given  by  H  =  <h  I  ♦  t  and 
x*(n».T»)  is  inverse  bremmstrahlung  absorption  coefficient.  Compression  effects  and  electron-ion 
energy  coupling  are  much  smaller  than  the  terms  included  in  (6.9).  and  are  ignored.  A  conservative 
flux-limited  formalism  is  used  for  the  heat  flux;  i.e.. 


t,  r  dTr  /dx  ■ 

f  e  ne  Te  vt 


(6.10) 


where  k ,  is  the  electron  thermal  conductivity  of  the  plasma  and  /,  is  the  flux  limiter,  which  is  set  to 
0. 1  here  in  order  to  ensure  that  Q,  does  not  exceed  its  classical  limit.  (For  the  simulations  presented 
here,  the  temperature  gradients  are  small  enough  that  the  heat  flux  remains  well  below  the  classical 
limit;  thus  the  magnitude  of  /,  has  little  effect  on  the  results.)  Equation  (6.9)  is  solved  by  a 
predictor-corrector  method  with  periodic  boundary  conditions.  This  can  cause  problems,  because  the 
net  energy  deposited  by  inverse  bremmstrahlung  cannot  be  lost  to  the  plasma  by  either  transverse  heat 
flow  losses  or  conversion  to  axial  plasma  kinetic  energy.  Therefore,  the  average  energy  gain 
H(z.t)  m  J dxHix.z .  t)/X,„  is  subtracted  from  the  source  term  of  (6.9)  in  order  to  keep  the  mean 
temperature  at  a  given  axial  distance  constant. 


For  steady-state  problems,  the  calculation  of  the  plasma  density  is  simplified.  The  momentum 
equation  for  constant  sources  ♦  :  and  T,  reduces  to  7u'1;nri  +  nf  aV  ♦  ;  =  o.  yielding  the  solution 


n  (.t .  c  ) 


cy(x  .z ) 


exp  -or  | 


3  ♦  :/dx 
c,*u  '  .z) 


(6.11) 


where  C'  is  a  constant.  The  temperature  T,  is  found  by  iterative  solution  of  Eq.  (6.9)  with  arbitrary 
time  steps,  performed  until  the  solution  converges  upon  itself.  These  separate  solutions  for  r,i.t.o 
and  ne(x,z)  are  iterated  alternately  to  obtain  a  convergent  steady-state  solution. 

B.  Results 

The  results  presented  here  compare  the  filamentation  behavior  exhibited  by  ISI  irradiation,  con¬ 
ventional  (quasi  near  field)  focusing.'*  and  the  RPS  technique.” 33  For  comparison  purposes,  we  start 
with  a  simple  well-defined  laser-plasma  model.  First,  the  plasma  is  chosen  to  be  homogeneous:  the 
constant  coupling  parameters  throughout  the  interaction  region  simplifies  the  analysis.  Second,  we 
ignore  absorption  in  the  plasma;  absorption  competes  with  filamentation.  and  can  mask  or  eradicate 
the  effect.  While  absorption  is  important  in  realistic  conditions,  it  unnecessarily  complicates  the 
interpretation  and  comparison  of  the  results.  Both  of  these  contraints  are  discarded  later  when  we 
look  at  a  "real  plasma". 

The  plasma  model  is  composed  of  a  homogeneous  250  tim  thick  CH  slab  at  nr„  =  0.5n  and 
Tt0  -  500eV,  irradiated  by  a  .25  jim  wavelength  laser  at  2  x  10u  W/cm3  Under  these  conditions, 
perturbation  theory38  predicts  that  both  the  thermal  and  ponderomotive  self-focusing  lengths  (for  the 
fastest  growing  modes)  are  comparable  (about  350X).  Thus,  the  interplay  between  these  two  mechan¬ 
isms  should  be  well  represented  here.  Note  that  this  plasma  model  is  selected  primarily  to  produce  a 
very  strong  filamentation  interaction;  while  useful  for  the  purpose  of  comparison,  it  gives  results  that 
are  more  pessimistic  than  those  of  actual  laboratory  plasmas,  as  we  will  see  later. 

Three  different  cases  are  investigated:  1S1  irradiance.  RPS  irradiance.  and  the  standard  quasi 
near-field  irradiance.  Of  these  three,  only  the  last  is  ambiguous  in  nature.  We  want  the  standard 
irradiance  case  to  represent  a  typical  high-power  laser.  In  practice,  the  structure  of  quasi  near-field 
illumination  can  vary  greatly,  depending  upon  the  laser  type,  the  power  level,  the  optical  engineering, 
the  attention  paid  to  beam  quality,  etc.  Here,  we  select  a  3:1  peak-to-average  as  characterizing  the 
typical  laser.  The  transverse  mode  spectrum  of  the  incident  plane  waves  is  chosen  to  be  fiat  for  all 
spatial  frequencies  between  T//?ipw  and  NSAK  =  k/2(F/number). 


Two  of  these  cases  under  scrutiny,  the  RPS  and  the  standard  laser  case,  require  only  the  steady 
state  analysis  because  the  intensity  patterns  remain  stationary  in  time.  First  consider  the  standard 
case,  where  the  peak  to  average  intensity  value  is  3:1  (the  standard  deviation  is  <j,  =  51).  Figure  14 
shows  the  intensity  profile  as  a  function  of  the  propagation  distance  into  the  plasma.  The  first  focal 
point  appears  at  z  =>  200  X,  and  additional  focal  points  appear  with  increasing  frequency  as  the  beam 
propagates  farther  into  the  plasma.  This  focal  distance  is  shorter  than  the  predictions  using  either  the 
ponderomotive  or  thermal  mechanisms  alone;  it  is  evidence  that  both  mechanisms  are  acting  together 
to  promote  filamentation.  The  filaments  that  are  created  appear  to  attract  one  another,  and  converge 
into  a  single  large  scale  mode  that  has  a  high  spatial  frequency  sub-structure.  This  behavior  is  caused 
by  the  interaction  of  the  thermal  and  ponderomotive  forces;  the  ponderomotive  force  is  primarily 
responsible  for  the  small  initial  filaments,  since  its  fastest  growing  mode  is  approximately  40X  in 
width.  Because  thermal  conduction  spreads  out  temperature  variations,  the  thermal  filamentation 
mechanism  only  acts  over  larger  distances.  The  fastest  growing  resolvable  thermal  mode  is  the  same 
size  as  the  simulation  dimension,  and  this  appears  to  be  responsible  for  the  large  scale  final  mode. 
Note  that  the  transverse  spatial  frequency  structure  of  the  irradiance  changes  fundamentally  as  the  fila¬ 
ments  are  formed;  after  the  first  few  filaments  form,  the  structure  is  spread  out  over  a  much  larger 


region  in  frequency  space.  At  first,  these  higher  spatial  frequency  structures  are  locally  coherent  and 
form  coherent  filaments.  Once  these  filaments  interact  with  one  another,  however,  the  intensity  pat 
tern  is  randomized  as  the  high  frequency  structure  becomes  incoherent  on  a  much  smaller  scale  The 
resulting  interference  pattern  seen  in  the  final  large  scale  mode  then  closely  resembles  the  initial 
structure  in  the  RPS  method.  Thus,  once  significant  filamentation  occurs  and  the  modes  interact 
(because  of  thermal  self- focusing),  further  filamentation  should  be  suppressed  in  accordance  with  RPS 
theory . " 

In  the  RPS  calculations,  there  is  a  moderate  amount  ot  suppression  or  filamentation  .1  >ne  uses 
small  F  number  i  <F  10)  optics.  Fig  15  illustrates  the  reduction  ot  filamentation  that  occurs  as  the 
F  number  decreases.  As  a  measure  of  the  degree  of  filamentation.  the  largest  intensity  <x.curnng  in 
the  plasma  is  used.  [Typical  maximum; average  values  ot  the  intensity  profile  incident  on  the  plasma 
are  as  follows:  (ii  for  the  standard  case  i/^  /,,  i  »„=  3.  mi  for  RPS  ‘I  w  /„  i.  =  4  to  b.  un  tor 
ISI  4  to  5  and  <  /  ^  /  .„  i- =  l  |  For  F  number  >F  lu.  the  maximum  intensities 

are.  on  average,  about  25  %  or  more  larger  than  for  F  number  <  F  10  There  may  be  better  suppres 
sion  below  f/2.5.  but  the  parabolic  wave  approximation  begins  to  break  down  there  and  the  simula¬ 
tion  results  are  not  trustworthy  Also,  the  trend  may  be  more  pronounced  at  higher  F  number  than 
shown  here.  The  runs  with  larger  F  number  do  not  include  the  longer  wavelength  modes,  which  are 
lumped  into  the  DC  term,  and  thus  the  peak  initial  intensities  are  somewhat  lower  than  migni  other 
wise  be  expected. 

The  ISI  simulation  was  run  for  350  coherence  times  ir  *  pseci.  and  most  of  the  results 
presented  here  were  averaged  over  this  entire  interval  The  most  important  result  is  the  time  aver 
aged  intensity  distribution  shown  in  Fig  16  This  distribution  exhibits  the  same  development  ot  high 
spatial  frequency  structure  found  tn  the  standard  laser  case,  but  without  the  development  ot  the  atten¬ 
dant  intensity  peaks.  The  suppression  of  these  peaks  appears  to  be  primarily  due  to  the  tact  that  the 
ISI  irradiance  structure  constantly  shifts  around,  not  allowing  any  single  filamentary  structure  to  per¬ 
sist  in  the  same  spot. 

However,  the  suppression  is  not  complete  on  the  shorter  time  scales  A  typical  instantaneous 
intensity  distribution  (Fig  1 73  shows  that  filamentary  structures  are  still  present,  although  at  a  more 
modest  level  than  in  the  standard  case.  To  compare.  Fig  18  plots  the  values  ot  /  ,u>  /  vs  propaga¬ 
tion  distance  for  the  standard  case,  along  with  the  maximum  value  of  the  time  averaged  ISI  intensitv 
i  /  ITUX).  and  the  time  average  of  the  maximum  ISI  intensity  <  /w  i  The  ISI  case  is  similar  to  the  con¬ 
trol  case,  except  for  two  significant  differences  First,  the  buildup  of  the  intensity  maximum  takes 
much  longer  with  ISI  (with  the  first  intensity  maximum  occurring  around  ;  =  500k.  xs  opposed  to 
;  *  200X  in  the  standard  case),  meaning  that  ISI  has  lengthened  the  self- focusing  distance  Second, 
the  peak  intensity  with  ISI  is  a  factor  of  two  lower  than  that  found  without  ISI  This  may  be  due  to 
the  smoothing  effect  in  the  time  average,  or  it  may  be  due  to  the  smaller  density  channels  made  in  the 
ISI  process.  (Note  that  although  the  averaged  peak  intensity  in  ISI  is  smaller  than  the  standard  case 
the  initial  averaged  peak  intensity  is  larger.) 

While  ISI  effectively  suppresses  filamentation  of  the  time-averaged  intensity,  the  behavior  of  the 
instantaneous  intensity  distribution  m  the  plasma  is  still  cause  for  concern,  as  many  plasma  instabili¬ 
ties  have  very  short  growth  times  (  -  psi,  and  can  therefore  react  to  the  peaks  seen  here  For¬ 
tunately.  laser  fusion  pellets  will  not  have  densities  -  0.5  or  sub-kilovolt  temperatures  extending 
over  hundreds  of  laser  wavelengths  As  an  example  of  a  more  realistic  plasma,  we  present  a  >imula- 
tion  taken  from  a  short-wavelength  pellet  interaction  used  in  an  earlier  study  *  The  target  is  a  DT 
pellet  driven  by  a  multi-megajoule  laser.  26  nsec  into  the  pulse,  after  a  significant  portion  ot  plasma 
volume  has  been  created.  The  temperature  and  density  profiles  used  for  this  example  are  plotted  in 
Fig  19  To  a  good  approximation  the  radial  convergence  can  be  ignored,  and  the  target  treated  as 
planar  The  major  differences  between  this  and  the  previous  model  are  the  inclusion  ot  absorption 
(approximately  50°£  of  the  light  is  absorbed  before  n  =  5 n  >.  higher  temperatures,  and  the  large 


volume  of  plasma  that  is  very  underdense  (and  therefore  more  weakly  coupled  to  the  light).  The 
results,  as  plotted  in  Fig.  20.  show  that  ISI  suppresses  filamentation  even  on  the  short  time  scales, 
and  negates  filamentation  as  a  problem.  On  the  other  hand,  the  results  with  the  standard  irradiation 
condition  (2.25:1  peak  to  average)  indicate  that  filamentation  will  be  a  nuisance  without  ISI  (Fig.  21 ): 
here,  a  filament  with  a  maximum  intensity  9  times  the  incident  average  intensity  has  been  formed. 

VII.  SUMMARY  .AND  DISCUSSION 

This  paper  has  described  the  ISI  technique  in  detail,  and  explored  its  potential  application  to 
direct-dnve  laser  fusion.  We  have  derived  the  ideal  ISI  profiles  in  the  complete  overlap  and  quadra¬ 
ture  configurations  [Eqs.  (2.17)  and  (2.23).  respectively),  and  examined  the  perturbing  factors,  such 
as  multibeam  interference,  laser  aberration,  and  plasma  filamentation.  which  can  introduce  nonumfor- 
miues  or  limit  the  controllability  of  those  profiles. 

Using  both  analytic  theory  and  numerical  simulations,  we  have  shown  that  ablation  pressure 
nonuniformities  due  to  the  rapidly  fluctuating  interference  structure  are  smoothed  by  two  independent 
mechanisms.  Thermal  diffusion  suppresses  the  higher  spatial  frequency  components,  in  accordance 
with  the  'cloudy  day"  model.  [E  g  see  Eqs.  (3.7)  or  1 3 .2 la.bi. |  If  this  mechanism  acted  alone  ias 
it  would  in  the  case  of  narrow  bandwidth  light)  the  standard  deviation  a„  would  be  approximately  pro¬ 
portional  to  iX.iiit/  D  i  [Expression  <3.25b>|.  which  favors  a  large  conduction  zone  width  a  and 
small  overall  F  number  r  D  Temporal  averaging  provides  an  additional  smoothing  factor 
J,  -  o. n‘  :  [Eq.  (3  19)),  which  favors  a  wide  laser  bandwidth  =  1/t  and  a  long  pulse  The 
standard  deviation  a,  =  i,  j,  resulting  from  the  combination  of  these  two  mechanisms  is  given 
approximately  by  expressions  (3  25a.b»;  this  result  shows  that  ablation  pressure  nonuniformities  - 
1  %  should  be  attainable  under  laser-plasma  conditions  of  interest  to  ICF  In  the  partial  1S1  scheme, 
there  is  a  stationary  component  in  the  interference  structure  due  to  the  coherent  repeating  steps 
Although  this  component  cannot  be  smoothed  by  temporal  averaging,  its  spatial  frequencies  are  suffi¬ 
ciently  high  that  they  are  almost  completely  eliminated  from  the  ablation  pressure  by  the  thermal 
smoothing  mechanism.  Thus,  partial  ISI  may  be  one  way  to  avoid  the  excessive  pulse  lengthening 
effects  that  would  otherwise  occur  in  large  aperture  systems  where  the  number  of  independent 
beamlets  could  exceed  St  If  necessary,  one  could  avoid  the  coherence  between  repealing  steps 
by  segmenting  or  multiplexing  the  laser  beam  itself,  driving  each  segment  by  either  separate  oscilla¬ 
tors  or  a  different  switched-out  portion  of  a  long  single  oscillator  pulse 

The  theory  has  shown  that  with  ISI  the  average  intensity  profile  at  the  target  remains  relatively 
insensitive  to  laser  beam  aberrations  whose  scale  lengths  are  larger  than  the  initial  beamlet  width 
D  In  practice,  this  means  that  an  incident  beam  aberrated  to  v„  nmes  its  diffraction  limit  would 
require  V,  >  j.V^  echelon  steps  to  ensure  adequate  control  over  the  target  profile.  The  aberration 
will  tend  to  misalign  and  broaden  this  profile  somewhat  [Eqs.  (4  10)  or  (4  12)|,  but  it  will  not  intro¬ 
duce  any  small-scale  structure  In  the  usual  case  where  >  >  i  and  the  aberration  arises  primarily 
from  random  phase  structure  on  the  beam,  the  relative  intensity  perturbations  due  to  misalignment  and 
broadening  are  R,  -  I  V,  [Eq.  (4  18)|  and  ,  .v,:  [Eq.  (4  20)|.  respectively  Most  of  the 

broadening  effect  is  controllable  because  it  depends  only  upon  spatial  averages  of  the  phase  structure 
over  the  enure  laser  aperture;  the  uncontrollable  component  is  A R:  ~  [Eq.  i4  22)1  Laser 

intensity  nonuniformities  can  also  contribute  to  perturbations  of  the  ISI  profile,  especially  in  the  qua¬ 
drature  configuration;  however  this  effect  also  depends  upon  a  spatial  average  over  the  aperture,  and 
is  expected  to  be  negligible  if  V,  >  >  I .  Our  results  show  that  it  should  be  possible  to  control  the  ISI 
intensity  perturbations  to  within  —  \%.  even  with  the  degree  of  aberration  expected  in  realistic  fusion 
drivers 

Two  dimensional  simulations  of  optical  propagation  using  the  SELFOCT  code  have  shown  that 
ISI  significantly  reduces  filamentation  in  the  underdense  plasma,  even  in  a  non-absorbing,  high  den¬ 
sity  (0.5  n  i.  low  temperature  (500  eVi  slab,  where  a  normal  laser  beam  would  experience  intense 
filamentation.  For  low  F  number  optics,  the  filamentation  also  appears  to  be  inhibited  by  the 


random-phase  shift  technique, ::  :j  but  to  a  lesser  degree  than  with  ISI.  In  simulations  at  0.25  ^m. 
including  inverse  bremsstrahlung,  higher  temperatures,  and  a  more  realistic  plasma  profile,  the  fila- 
mentauon  is  virtually  eliminated  with  ISI.  (Fig.  20) 


An  important  concern  tor  ISI  is  the  question  of  how  broad  bandwidths  or  transient  nonuniformi¬ 
ties  would  affect  certain  laser-plasma  instabilities  (e.g.  parametric  instabilities),  which  can  respond 
much  faster  than  the  hydrodynamics.  Earlier  experiments  at  1.054  indicated  that  although  broad 
laser  bandwidth  suppressed  stimulated  Brillouin  scattering  in  the  underdense  plasma,  as  predicted 
theoretically,  the  bandwidth  appeared  to  enhance  a  critical  region  instability  that  produces  fast  elec¬ 
trons.’940  In  the  NRL  experiments  at  1  micron.10  this  effect  was  observed  both  with  and  without  the 
ISI  echelons.  The  bandwidths  used  in  these  experiments  were  comparable  to  the  expected  ion- 
acoustic  frequency  for  the  parametric  decay  instability,  and  the  associated  temporal  mode  beating  may 
have  been  seeding  this  critical  region  instability.  More  recently, :o  the  experiments  were  repeated  at 
0.527  nm  to  determine  whether  the  use  of  a  shorter  laser  wavelength,  which  is  strongly  absorbed 
before  it  reaches  the  critical  region,  would  inhibit  the  fast  electron  generation.  This  indeed  occurred: 
in  fact,  the  experiments  showed  that  with  the  combination  of  ISI  and  shorter  wavelength,  the 
enhanced  hot  electron  generation  observed  earlier  was  replaced  by  a  virtual  quenching  of  hot  electron 
production.  In  summary,  the  NRL  experiments  are  consistent  with  the  following  model:  (1)  The 
smooth  irradiance  profile  obtained  with  ISI  echelons  and  broad  bandwidth  suppresses  the  effects  of 
hot  spots,  preventing  for  example  the  seeding  of  self  focusing.  (2)  The  combination  of  this  smoothing 
with  the  stabilizing  effect  of  the  bandwidth  itself  suppresses  underdense  parametric  instabilities.  (3) 
The  shorter  laser  wavelength  prevents  light  from  reaching  the  critical  region,  and  thereby  suppresses 
broadband  enhancement  of  hot  electrons. 

In  the  ISI  and  other  related  schemes  discussed  so  far.  the  beam  is  divided  up  by  optical  arrays 
placed  at  the  output  of  the  laser,  where  coating  damage  and  optical  complexity18  ’6  may  be  possible 
issues.  This  restriction  on  the  location  of  the  ISI  arrays  stems  from  self  focusing  effects  in  glass 
lasers,  and  from  the  necessity  of  maintaining  spatial  coherence  in  any  harmonic  conversion  crystals. 
Near-fietd  nonuniformides  associated  with  an  array  of  beamiets  would  seed  self  focusing  if  one 
attempted  to  amplify  those  beamiets  in  a  multistage  glass  laser.41  In  KrF  lasers,  this  restriction  is  no 
longer  necessary  because  the  amplifying  medium  is  gaseous,  intensities  are  low  (typically  <  10 
MW/cnr).  and  harmonic  conversion  is  not  required.  The  spatial  incoherence  could  then  be  induced  by 
echelons42  or  some  alternative  technique  at  a  low  energy  stage  within  the  laser,  and  optically-relayed 
through  the  remaining  stages  of  the  amplifier  chain.  For  example,  one  could  use  reflection  from  a 
so-called  plasma  “spatial  filter''  to  generate  the  spatially- incoherent  light.43  Although  this  technique  is 
capable  of  producing  a  smooth  focal  spot,  43  44  it  is  not  clear  whether  there  would  be  sufficient  control 
over  the  spatial  coherence  properties,  or  sufficiently  rapid  temporal  averaging  of  the  interference 
structure  to  ensure  that  the  target  profile  is  controllable.  Alternatively,  one  could  use  the  scheme 
shown  in  Fig.  22,  in  which  a  beam  of  statistically-homogeneous  incoherent  light  is  generated  by 
broadband  amplified  spontaneous  emission  (ASE)  in  a  mirrorless  or  spatially  multimode  oscillator.21 
This  beam  traverses  a  variable-density  absorber,  whose  spatial  transmission  function  Fix)  is  the 
desired  target  beam  profile.  It  is  then  focused  onto  the  entrance  pupil  of  an  optically-relayed  KrF 
laser  chain.  Because  the  entrance  pupil  lies  at  the  Fourier  transform  plane  of  the  absorber,  the  optical 
information  needed  to  reproduce  F(x)  at  the  target  will  be  transported  through  the  laser  by  a  multitude 
of  small  coherence  zones,  rather  than  any  large  whole  beam  structure:  thus  the  coherence  zones  play 
essentially  the  same  role  as  the  beamiets  do  in  the  conventional  ISI.  The  target  beam  profile  will  be 
insensitive  to  laser  imperfections  as  long  as  the  coherence  zones  remain  small  in  comparison  to  the 
transverse  scale  lengths  of  the  phase  aberration  or  amplifier  gain  nonuniformities.  It  will  be  insensi¬ 
tive  to  amplifier  saturation  if  the  coherence  time  tc  remains  short  in  comparison  to  the  KrF  relaxation 
time  (>l/u).  Theoretical  analysis,  along  with  a  preliminary  experimental  test  using  a  small  KrF 
oscillator-amplifier  setup,  has  demonstrated  that  this  technique  is  indeed  capable  of  producing  a 
smooth,  controllable  focal  profile,  even  under  non-optimum  conditions.21  Additional  results  on  a 
larger  KrF  system  will  be  reported  in  later  publications. 
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APPENDIX 


This  Appendix  calculates  the  moments  of  the  target  profile  that  the  laser  would  produce  without 
the  ISI  echelons.  The  coherent  (but  aberrated)  amplitude  in  this  plane  is  given  by  the  diffraction 
integral20 
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Substituting  expression  (2.1)  for  EL{x‘.n,  one  obtains  the  intensity 
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where  Cc  -  b2c/%*.  From  Eq.  (2.5)  and  the  definition  b2  =  &*/cD; ,  one  obtains 
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The  centroid  of  Ic{x.z)  is  defined  by  Eq.  (4.7) 
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Performing  the  integration  over  x,  which  gives  3(x’  -x  ’),  then  integrating  by  parts  over  x  .  one 
finally  obtains 
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which  reduces  to  (4.6a)  for  :  <  </ .  The  second  spatial  moment  is  defined  by 
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Following  steps  similar  to  those  used  on  (A4).  one  obtains  expression  (4.6b)  for  :  <  <  / 
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FOCAL  IRRADIANCE 

Figure  I  —  Schematic  of  the  ISI  concept  in  two  dimensions.  The  broadband  laser  beam  of  coherence  time  /,  is 
broken  up  into  mutually  incoherent  heamlels  it  the  echelon-induced  time  delay  increments  satisfy  ,  —  t„  >  i:  . 

The  lens  combines  these  bcamlcls  at  its  local  point  to  give  a  profile  that  has  a  smooth  lime  average  for  times  > 
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UASI  NEAR  FIELD  FOCUSING  C 
ABERRATED  LASER  BEAM 


FLUENCE  AT  THE  ABSORPTION  REGION  ABLATION  PRESSURE,  ASSUMING 

THERMAL  SMOOTHING  WITHd=  50/im. 

Figure  8  Numerical  simulation  ol  (lie  conventional  tocusing  geometry,  where  the  target  is  placed  in  the  quasi 
near  held  ol  the  5  X  1)1  alienated  lasci  beam.  The  ablation  pressure  retains  large  nonunilorniitics  in  spile  ol 
the  thermal  smoothing 


I.S.I.  WITH  BROADBAND  LIGHT 


FOCAL  POINT  FLUENCE  WITH  THE 
ECHELONS  IN  PLACE,  ANO  A  BROADBAND 
LASER  INTEGRATED  OVER  TIME  INTERVAL 

r  -  100  t,.. 


ABLATION  PRESSURE,  ASSUMING 
THERMAL  SMOOTHING  WITH  d=  50^m. 


Figure  10  —  Numerical  simulation  of  ISI  with  broadband  527  nm  light,  showing  the  quadrature  intensity  and 
ablation  pressure  profiles  averaged  over  100  coherence  times.  The  peak/vailey  deviation  from  the  ensemble- 
averaged  ablation  pressure  is  now  2.4%. 
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Figure  15  —  Comparison  of  maximum  intensity  achieved  inside  plasma  as  a  function  of  incident  lens 
/  number,  for  a  variety  of  input  intensity  profiles. 


Figure  18  —  Comparison  of  peak  values  of  ISI  and  standard  laser  irradiation  vs  propagation  distance  in  plasma. 
The  solid  line  is  the  “typical”  laser  (3:1  peak-to-average  incident  intensity  variations),  the  dotted  line  is  the 
time  averaged  ISI  intensity  peaks  ({/„»«/),  and  the  dashed  line  is  the  maximum  of  the  time  averaged  ISI 
distribution  (( l) 
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Hgute  22  Schematic  diagram  ol  a  spuiially  incoherent  broadband  laser  system  lor  smooth,  controllable 
illumination  ot  fusion  targets  I  he  line  line  profiles  show  typical  instantaneous  intensities  |  £  | ’,  while  the 
heavy  lines  show  the  lime  avciagcd  intensities  |  £  |  * 
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